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I.  INTRODUCTION 


A.  BACKGROUND 

Solid  propellant  rocket  motor  design  is  a  very  complex 
process  which  often  requires  many  iterations  in  order  to 
accomplish  a  successful  design.  Usually,  the  first  step  in 
the  design  process  it  to  obtain  accurate  input  data  on  mis¬ 
sile  system  restraints.  Examples  of  these  restraints  are 
ballistic  requirements  (which  can  include  thrust,  burning 
time,  total  impulse)  and  mission/vehicle  requirements  (which 
can  include  motor  length,  diameter,  allowed  hazard  and  plume 
properties,  altitude,  etc.).  These  parameters  are  often  re¬ 
ferred  to  as  independent  parameters.  Through  the  use  of 
these  independent  parameters,  the  required  propellant  burn 
time  and  properties  can  be  roughly  estimated.  Then,  theo¬ 
retical  motor  performance  and  losses  can  be  estimated.  The 
thrust-time  requirements  together  with  the  expected  struc¬ 
tural  loads  permit  initial  selection  of  a  grain  design. 

The  next  step  of  solid  propellant  rocket  motor  design  is 
to  use  the  prooellant  properties  and  grain  design  tu  predict 
the  pressure-  and  thrust-time  profiles  and  the  total  impulse 
of  the  motor.  Propellant  properties  and  grain  design  are 
then  iterated  until  these  independent  variable  requirements 
are  met,  together  with  acceptable  propellant  stress  levels. 
A  final  step  in  the  missile  design  process  is  to  use  the 


thrust-time  nrofile  in  an  external  aerodynamics  code  to 
ensure  that  range  and  terminal  velocity  requirements  are 
met . 

During  the  past  decade,  there  has  been  significant  pro¬ 
gress  in  the  automation  of  solid  rocket  motor  design.  It  is 
mainly  because  of  the  extensive  use  of  digital  computers, 
which  greatly  facilitate  the  iterative  type  computations 
mentioned  above  (Ref.  1).  In  recent  years,  many  of  the 
design  processes  have  been  accomplished  using  separate  com¬ 
puter  codes,  i.e.,  adiabatic  equilibrium  combustion  codes, 
codes  which  estimate  performance  losses,  grain  geometry 
optimization  codes,  etc.  These  codes  have,  in  some  cases, 
been  combined  into  much  larger  design  codes,  and  they  are 
often  difficult  and  expensive  to  run.  Recently,  PC  versions 
of  some  rocket  motor  design  codes  have  become  available. 
Hov^/ever,  they  must  still  be  used  individually,  and  were  not 
all  available  for  use  on  the  sane  computer.  The  use  of 
these  codes  was  also  very  inconvenient  in  that  the  same 
parameters  often  had  different  variable  names  in  the  various 
codes.  The  user  had  to  transfer  or  convert  the  data  from 
one  code  to  the  next. 

In  addition,  a  recent  AGARD  publication  has  provided 
recommended  emirical  equations  which  can  be  used  for  rapid 
and  accurate  estimation  of  performance  losses  for  aluminized 


propellants . 


WIlW'ALWi.'*] 
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B.  OBJECTIVES 

The  first  objective  was  to  outline  a  suitable  "Solid 
Propellant  Rocket  Motor  Preliminary  Design  Procedure"  which, 
of  course,  is  not  a  unique  procedure.  The  second  objective 
was  to  create,  modify,  and/or  combine  the  PC  versions  of  the 
necessary  design  codes  into  one  interactive  code  which  could 
be  conveniently  used  for  preliminary  aesign  on  an  IBM  PC/AT 
microcomputer.  The  final  objective  was  to  demonstrate  the 
program  by  designing  a  tactical  solid  propellant  rocket 
motor  combustor. 

C.  APPROACH 

The  approach  to  achieve  the  objectives  was  first  to 
detail  an  acceptable  preliminary  design  process.  The  algo¬ 
rithm  of  the  design  process  was  then  followed,  utilizing 
several  existing  codes  and  loss  estimation  methods  as 
appropriate,  to  obtain  an  integrated  design  code  for  use  on 
the  IBM  PC/AT.  The  programs  to  be  used  were  (1)  "MICROPEP" 
(a  microcomputer  version  of  a  propellant  evaluation  program 
from  China  Lake,  Naval  Weaoons  Center),  which  calculates  the 
adiabatic,  equilibrium,  combustion  process,  and  theoretical 
motor  performance;  (2)  "GRAINS"  (a  radial  spoked  star  grain 
design  program  from  CSD);  (3)  "ROCKET  (a  Lockheed  internal 
ballistics  program),  which  is  a  one-dimensional  flow  code 
which  predicts  pressure  and  thrust  versus  time;  and  (4) 
"FLYIT"  (a  flight  simulation  program  from  China  Lake,  Naval 
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Weaoons  Center),  which  perforns  simplified  trajectory  calcu¬ 
lations.  In  addition,  a  code  based  upon  AGARD-a,r-230  was 
written  for  estimation  of  performance  losses  for  aluminized 
propellants.  Capabilities  of  the  final  program  were  then  to 
be  demonstrated  by  considering  the  design  of  a  tactical 


motor . 


II.  SOLID  PROPELLANT  ROCKET  MOTOR  DESIGN  PROCESS 


A.  GENERAL 

The  rnethod  of  solid  rocket  motor  design  adopted  in  this 
work  is  shown  in  Figure  1.  It  was  assumed  that  the  required 
information  on  fundamental  propellant  properties  (burning 
rates,  temperature  sensitivity,  etc.)  was  known.  In 
addition,  the  design  process  discussf^d  herein  did  not 
include  any  stress  analysis  for  the  propellant  grain  or 
case . 


B.  INPUTS 

The  inputs  that  are  required  to  satisfy  mission  objec¬ 
tives  (the  independent  parameters)  were  as  follows: 

1 .  Ballistics  Performance  Requirements 

a.  Average  thrust  (F)  or  thrust-time  profiles, 

b.  Burn  time  (t^),  or 

c.  Total  impulse  (It^* 

These  parameters  are  interrelated  as  shown  in 


Figure  2.  Thus  =  j  Fdt 

o 


F  =  It/tb 

2.  Mission/Vehicle  Constraints 


2.1) 


(2.2) 


These  requirements  are  directly  related  to  the  motor 


design  and  include: 


Si 


a.  Envelope 


The  allowable  envelope  (motor  length  L,  motor 
inside  diameter  D)  is  the  fundamental  constraint  on  the 
motor  geometry,  which  can  be  ,  used  to  estimate  the  reguired 
burning  rate  (r)  and  total  volum.e  available  (V^)  for  the 
propellant  by  following  relations; 

r  =  (D/4)/tb  (2.3) 

(i.e.,  a  web  fraction  of  approximately  .5) 

Va  =  (  ttDV4)  L  (2.4) 

b.  Maximum  Expected  Operating  Pressure  (MEOP) 

The  MEOP  is  set  by  the  structural  limit  of  the 
motor  case.  The  maximum  operating  pressure  at  maximum  pro¬ 
pellant  temperature  (Pmaxtp^  less  than  the  MEOP  by  the 

factor  of  reproducibility  tolerance  (see  Figure  2).  Pmaxtp 

can  be  initially  estimated  from  the  following  equation 

Pmaxtp" (1-reproducibility  tolerance)  (2.5) 

c.  Temperature  Range  and  Altitude 

The  operating  propellant  temperature  range 

effects  the  selection  of  the  propellant  through  the  required 
burning  rate  temperature  sensitivity.  The  design  altitude, 
together  with  the  nomiiial  chamber  pressure  determine  the 


exhaust  nozzle  area  ratio. 


j« *. ji*. ji^  W>  •  ••  k  \  . 
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d.  Allowed  Hazard  Rating  and  Allowed  Exhaust 
S ignature 

The  hazard  rating  requirements  have  to  do  with 
whether  the  propellant  is  detonable  or  non-detonable . 
Modern  high  energy  propellants  are  often  detonable.  The 
allowed  exhaust  signature  affects  the  propellant  selection 
through  requirements  on  exhaust  smoke  and  radiation. 


C.  PROPELLANT  SELECTION 

From  the  ballistics  performance  and  mission/vehicle  con¬ 
straints,  it  is  possible  to  make  a  preliminary  selection  of 
the  propellant  type  that  is  most  likely  to  give  the  reguired 
performance,  internal  ballistics,  flame  temperature  and 
mechanical  properties,  as  well  as  the  necessary  storage 
stability,  the  best  hazard  properties  and  exhaust  signature 
{Ref.  2).  A  particular  propellant  will  have  known  (or 
estimated)  burning  rate  (r),  temperature  sensitivity 
density  (pp)f  and  composition.  Then  the  nominal  chamber 

pressure  (P^  )  can  be  calculated  by  using  the  equation 

nom 

(Ref.  3): 


^c  "  Pmaxtp  ("^nom  Tj^ax^) 

nom  ^ 


(2.6) 


D.  THEORETICAL  PERFORMANCE  COMPUTATION 

Following  the  propellant  selection,  theoretical  per¬ 
formance  computations  can  be  made  for  the  specified  chamber 
and  ambient  pressures.  A  program  such  as  PEPCODE  (Ref.  4), 


or  the  PC  version  "fllCROPEP",  will  provide  data  such  as 
thrust  coefficient,  characteristic  exhaust  velocity,  shift¬ 
ing  and  frozen  equilibrium  performance  results,  optimum  noz¬ 
zle  expansion  ratio,  mole  fraction  of  condensed  species  in 
the  exhaust,  etc. 

E.  NOZZLE  GEOMETRY  AND  MATERIAL  SELECTION 

The  information  on  theoretical  combustion  temperature, 
pressure,  gas  composition,  and  burn  time,  together  with 
geometry  constraints,  permits  the  nozzle  configuration  and 
material  to  be  selected. 

F.  PERFORMANCE  LOSS  ESTIMATION 

It  was  assumed  that  the  delivered  specific  impulse  could 
be  written: 

Isp  ^sPth 

where  1=,^  is  the  theoretical  specific  impulse  calculated 
®Pth 

by  "MICROPEP";  and  nc^  ^C*  thrust  coefficient 

efficiency  and  characteristic  velocity  efficiency  (or  com¬ 
bustion  efficiency),  respectively. 

nr  was  computed  by  summing  the  effects  from  the 
^f 

following  losses 

1.  Divergence  losses  ( £dIV^ ' 

2.  Two-phase  flow  losses  (e^p), 

3.  Boundary  layer  loss  (epL^' 
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5.  Submergence  loss  (esuB^' 

6.  Erosion  loss  (geROS^  ■ 

Thus,  nc^  =  (  ^DIV*'^TP‘^^BL'*'^KIN'''^SUB‘''^ER0S  (2.8) 

The  combusion  efficiency  depends  primarily  on  the 
residence  time  (tj-),  and  the  latter  can  be  expressed  as 

t^  =  (Vol/m)  M/12  R  T^)  (2.9) 

where 


- 

Vol 

is 

the 

instantaneous  volume  of  the 

combusor  ( in ^ ) , 

- 

Pc 

is 

the 

chamber  pressure  (Ibf/in^), 

- 

Tc 

is 

the 

chamber  temperature  ( °R) , 

- 

M 

is 

the 

molecular  weight. 

- 

R 

is 

the 

universal  gas  constant  (ft 

lbf/lbn,oleOR)  ,and 

- 

• 

m 

is 

the 

mass  discharge  rate  (Ibm/sec). 

For 

example , 

Figure  3  (Ref.  2)  is  a  plot  of 

specific  impulse 

versus  average  residence  time  for  aluminized  propellants. 
It  shows  that  the  residence  time  must  be  somewhere  between 
10  to  15  msec  to  obtain  a  combustion  efficiency  between  90% 
and  95%.  Increasing  the  free  volume  of  the  combustion 
chamber  will  improve  the  residence  time.  In  this  work  the 
initial  volume  was  used  in  the  residence  time  calculation, 
i.e.,  the  worst  case. 

G.  THROAT  AREA  AND  EXIT  AREA 


Accurately  sizing  the  throat  and  exit  areas  is  one  of 
the  major  goals  of  the  solid  propellant  rocket  motor  design. 


0  -  5  10  15  20  25  30  35  40  45  50 

Average  residence  tiiTie,  nsec 


FIGURE  3. 

EFFECT  OF  RESIDENCE  TIME  ON  ^ sp  EFFICIENCY (Ref 
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These  areas  are  dependent  on  several  other  variables. 


such  as  thrust  coefficient  efficiency,  combustion 
efficiency,  specific  impulse,  and  etc..  In  the  sizing 
process,  it  is  convenient  to  start  by  using  the  theoretical 
values.  The  effects  of  losses  can  then  be  handled  itera¬ 
tively  until  the  solution  converges  to  the  final  areas.  The 
rest  of  this  section  will  present  the  basic  algorithm  for 
the  throat  and  exit  area  sizing  process. 

The  basic  definition  of  specific  impulse  can  be  ex¬ 
pressed  as  (Ref.  3): 

Isp  =  F/m  (2.10) 


where 


Then 


Ath 


F _ 

Cfth  nCf  Pc 


(2.11) 


(2.12  ) 


The  theoretical  throat  area  can  be  estimated  by  using 

Cf  for  shifting  ecuilibrium  flow  (from  ’MICROPEP"  -  the 
th 

theoretical  performance  computation).  The  losses  which  are 
calculated  to  obtain  depend  upon  the  output  from 

"MICROPEP"  and  the  throat  area.  Thus,  iteration  is  reguired 
to  obtain  both  the  losses  and  throat  area. 

In  summary,  the  calculation  can  be  conducted  in  the 
following  steps: 
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1.  Calculate  using  equation  (2.11). 

2.  Calculate  ng  discussed  above. 

3.  Calculate  A^h  using  equation  (2.12). 

4.  Repeat  steps  2.  and  3.  until  the  values  for  throat 
area  and  nr  converge. 

In  the  current  program  the  iteration  was  terminated  when 
there  was  less  than  a  1%  change  in  A^-h  iteration. 


The  exit  area  can  then  be  calculated  from 


Ae  ~  ^th  ^ 


(2.13) 


where  e  is  the  exit  to  throat  area  ratio.  e  is  calculated 
in  "niCROPEP",  assuming  that  the  nozzle  exit  pressure  is 
equal  to  the  local  ambient  pressure  (one-dimensional,  ideal 
flow  is  assumed ) . 


H.  GRAIN  DESIGN 

The  primary  purpose  of  grain  design  is  to  shape  the  pro¬ 
pellant  grain  such  that  the  rocl<:et  motor  can  deliver  the 
required  performance  (thrust-time  or  pressure- t ime  profile). 

In  order  to  be  able  to  select  the  general  conf igurat ion 
of  the  grain  geometry,  several  parameters  must  first  be 
determined;  the  web  fraction  {V!^)  ,  volumetric  loading  frac¬ 
tion  (Vj_),  and  the  length-to-diameter  ratio  (L/D)  (Ref.  5). 
L/D  is  known  from  the  input  data.  Web  fraction  is  the  frac¬ 
tion  of  the  motor  radius  that  is  filled  with  propellant. 


Thus , 


Wf  - 


(  D/2) 


2.14  ) 


The  propellant  burning  rate  (r)  is  generally  expressed  in 
the  form; 

r  =  apc  (2.15) 
where  a  and  n  are  empirical  constants.  They  are  assumed  to 
be  known  for  the  selected  propellant.  The  volumetric  load¬ 
ing  fraction  is  defined  as  the  ratio  of  the  volume  occupied 
by  the  propellant  (Vj^)  to  the  total  volume  available  in  the 
motor  for  propellant  (Vg).  ThuS/ 


where 


Va  = 


(2.16) 


(2.17) 


The  required  propellant  volume  is  found  as  shown  below 


Wp 

it  “  ^^b  “  isp'^^b  “  isp  *^b  “  ispWp 


(2.18) 


Thus,  Wp  =  It/Isp  (2.19) 

This  yields  the  required  propellant  weight  (Wp)  and 

Vn  =  (2.20) 

Pp 

where  pp  is  the  density  or  specific  weight  of  the  selected 
propellant.  With  Wf ,  VI,  and  L/D  known,  an  initial  selec¬ 
tion  of  grain  qeometry  can  be  made  (for  example,  see  Ref. 
5).  For  a  selected  grain  configuration,  the  web  versus 


» 


^  * 


burning  surface 


area  (which  is  the  data 


needed  for 


thrust  and  pressure  versus  time  calculations)  can  be 
computed . 


I.  THRUST  AMD  PRESSURE  VERSUS  TIME  CALCULATIONS 

The  thrust  and  pressure  versus  time  calculations  are 
made  using  conservation  of  mass  for  the  combustor ( Ref .  3). 
The  mass  addition  from  the  burning  propellant  (nig)  is 
equal  to  the  sum  of  the  mass  accumulated  in  the  combustor 
and  the  mass  exhausted  from  the  nozzle  (m^j)  as  shown  in 
Figure  4.  Thus, 


mq  =  dM/dt  +  m^, 


(2.21) 


mQ  is  the  mass  of  gas  generated  from  the  burning  propellant. 


™g  "  I  Pp-^b*^  “  IPp^b^i^c 


(2.22) 


The  summation  is  over  all  propellant  grains  that  are  burning 
at  the  same  time. 


is  the  mass  accumulating  in  the  combustor. 


dfl  d  ,  „  ,  \  avoi  ,  ,,  ,  ^  ^^g 


dVol 


where 


then 


dVol  v’.  ~  ^ 

—  =  IrAb  and  _  =  _  ^ 


Pq  dVol  vol  dPc 


^2.23) 


(2.24 


/ 

/ 

7  /  / 

□M/dt 

/ 

/ 

/  / 

is  the  mass  exhausted  from  the  nozzle.  For  chocked  flow. 


m. 


rhn 


"^th  9c 


Combining  equations  2.21  through  2.25  results  in 


1 Pp^b^^c 


,  ^c^  dVolj_,  Vol  x'^^C  ^  r^C-^th^Ci 


RT, 


Solving  for  dP^/dt 


dP, 


dt 


=  ( 


Vol 


[rTc(I  Pp^b^^c 


^c'^th^c  X 
- rTf - > 


p  , dVol 
dt 


(2.25) 


(2.26 


)]  (2.27) 


Equation  2.27  can  be  numerically  inteqrated  to  give  P^^  (  t )  . 
The  thrust-time  behavior  can  then  be  obtained  from 

(2.28) 


F  (t)  =  nc^  Cf ( t)Pc( t)Ath( t) 


J.  TERMINAL  VELOCITY  AND  RANGE  CALCULATIONS 

The  last  step  of  solid  propellant  rocket  motor  design  is 
to  incorporate  the  designed  motor  into  the  desired  airframe 
and  make  a  flight  simulation.  The  purpose  of  this  flight 
simulation  is  to  verify  that  the  motor  can  produce  the 
desired  range  and  terminal  velocity. 


2':’ 


.Aiiairfe 


».  -w  -.  ■%.  "w  .  % 


III .  MI CROP EP 

"MICROPEP"  is  a  microcomputer  version  of  the  propellant 
Evaluation  Program  (PEP),  published  in  Naval  Weapons  Center 
Report  NVJC-TP-6037  ,  "Theoretical  Computations  of  Eguilibrium 
Compositions,  Thermodynamic  Properties,  and  Performance 
Characteristics  of  Propellant  Systems,"  by  D.  R.  Cruise, 
April  1979  (Ref.  4).  The  code  (MICROPEP)  can  handle  12  dif¬ 
ferent  input  ingredients  and  200  combustion  species.  It 
performs  the  calculations  for  the  thermodynamic  properties 
and  oerformance  of  propellant  systems. 

A.  BASIC  ASSUMPTIONS 

The  basic  assumptions  of  "MICROPEP"  are  as  follows: 

1.  Equilibrium  adiabatic  combustion  at  the  nozzle 

entrance . 

2.  One-dimensional,  isentropic  expansion  in  the  nozzle  in 
which  the  flow  can  be  either  a  shifting  equilibrium 
flow  (assuming  chemical  equilibrium  is  maintained 
throughout  the  expansion  process)  or  frozen  composi¬ 
tion  flow  (assuming  the  chemical  composition  of  the 
the  flow  remains  the  same  throughout  the  expansion 
process ) . 

3.  Gases  were  considered  to  be  perfect. 

4.  Nozzle  exit  pressure  is  assumed  to  be  equal  to  the 
ambient  pressure. 
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INPUTS 


Inputs  required  by  "MICROPEP'  are: 

Composition  and  mass  of  each  ingredient.  Normally  the 
weights  of  ingredients  are  chosen  to  add  up  to  100  q. 
NOTE:  Ingredients  may  be  input  by  the  user  or  they 

may  be  selected  from  an  existing  data  file.  In 
the  latter  case  ingredient  serial  numbers  are 
used . 

For  assistance,  use  README,  which  should  be  on 
the  PEPCODED.DAF  Data  File. 

Combustion  pressure. 

Nozzle  exit  pressure.  This  defines  the  limits  for  the 
expansion  process. 

OUTPUTS 

The  program  calculates: 

Combustion  product  composition,  properties  ,  and 
temperature  at  the  nozzle  entrance. 

Nozzle  exit  product  composition,  properties,  and 
temperature  for  two  cases: 

a.  Shifting  equilibrium  flow 

b.  Frozen  composition  flow. 

Required  throat  area  and  exit  area  with  no  losses,  for 
the  same  two  cases. 

Theoretical  specific  impulse,  thrust  coefficient,  and 
characteristic  exhaust  velocity  for  both  shifting 
equilibrium  and  frozen  composition  flow. 


i 


IV.  LOSS  MECHANISMS 


For  estimation  of  the  performance  losses,  a  set  of 
empirical  formulas  suggested  by  the  Propulsion  and 
Energetics  Panel  of  AGARD  in  AGARD-AR-230  (Ref.  6)  were 
used.  Similar  methods  can  also  be  found  in  NASA  Report  SP- 
8039  (Ref.  7)  "Solid  Rocket  Motor  Performance  Analysis  and 
Prediction."  It  is  assumed,  as  in  eguation  2.7,  that  the 
combustion  efficiency  (hc*^  thrust  coefficient  effici¬ 
ency  (  Tic  ^  treated  ^  independently.  The  empirical 
formulas  are  valid  only  for  aluminized  propellants. 


A.  THRUST  COEFFICIENT  EFFICIENCY 

The  thrust  coefficient  efficiency  (hc^)  can  be  expressed 
as  in  eguation  2.8 


ric^  =  1- (  eDIV'''^TP‘^^BL'*’^KIN''’CSUB'''CEROS 

Empirical  equations  for  each  of  the  losses  (e's)  will  be 
presented  and  briefly  discussed  below.  More  details  can  be 
found  in  Reference  6. 

1 .  Divergence  Loss  (^DIV) 

This  is  the  loss  due  to  the  radial  velocity  com¬ 
ponent  of  the  gas  at  the  nozzle  exit  plane 


GDIV  "  [1-Cos  (- 


(4.1) 
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a  IS  the  nozzle  half  angle  and  9ex  is  the  exit  angle  for  a 
contoured  nozzle.  For  a  contoured  nozzle,  a  is  measured  to 
a  line  drawn  from  the  exit  and  tangent  to  the  wall  at  the 
throat . 

2 .  Kinetics  Loss  (‘•KIN) 

This  is  the  loss  (or  reduction  in  performance  from 
the  value  of  shifting  equilibrium  flow)  associated  with  not 
attaining  chemical  equilibrium  throughout  the  nozzle  expan¬ 
sion  process.  It  is  estimated  to  be  1/3  of  the  fractional 
difference  between  the  specific  impulse  for  shifting  equi¬ 
librium  (Ion  )  snd  frozen  composition  (Ion  )•  Thus, 


£KIN 


=  33.3  [1  -  ] 

^sPs 


(4.2) 


Both  ^spg  and  ^spf  are  obtained  from  'MICROPEP." 

3 .  Boundary  Layer  Loss  ( ^BL) 

The  boundary  layer  loss  (egL^  expressed  as 


p  0.8 

,  f_c _ 

'1  ^  0.2 
Dth 


(1+2  exp(-C2 


tp)]  [l  +  0.016(e-9 


(4.3) 


This  loss  accounts  for  both  reduced  flow  areas  and  transient 
heat  loss.  In  equation  4.3,  P^,  is  in  psi,  Dj-p  is  in  inches. 
The  Cl  and  Co  coefficients  are 


for 

ordinary  nozzles 

^1  ^ 

0. 

.00365 

C2  = 

0. 

.000937 

for 

steel  nozzles 

Cl  = 

0. 

.00506 

C2  = 

0. 

.0 

4 .  Two-Phase  Flow 

Loss  ( ^TP ) 

This  loss  is  due  primarily  to  the  velocity  lag  of 
the  condensed  species  of  the  combustion  products  as  the 
mixture  passes  through  the  exhaust  nozzle. 


where 


STP 


=  C- 


Mf  C4 


p^O.15  e0*08 


(4.4) 


•^th 


a .  D 


th 


is  the  diameter  of  throat  in  inches 


b.  Dp  is  the  mean  AI2O3  particle  diameter  in 
microns  and  is  calculated  from 


Dn  =  3.39 


c . 


,p  - 

is  chamber  pressure  in  psia 


e.  Mf  is  the  mole  fraction  of  consdensed  phase  in 
moles/100  gm  of  reactants. 
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C's  are  dependent  on  the  diameter 

C4  =  0.5 

of 

the 

throat 

Dth  <1: 

C3  =  9  C5 

= 

1 

C6 

=  1 

1  <  Dth  < 

2;  C3  =  9  C5 

= 

1 

C6 

=  0.08 

Dp^>2  and 

and 

Dp<4 ;  €3=13.4  C5 

0.8 

C6 

=  0.8 

• 

j-  * 

4<Dp<8 : 

C3  =  10.2  C5 

0.8 

C6 

=  0.4 

v‘3 

«  ■  t 

Dp}^>2  and 

Dp>8;  C3  =  7.58  C5 

= 

0.8 

C6 

=  0.33 

^y:\ 

• 

■ 
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Submergence  Loss  ( ^SUB ) 


This  loss  is  expressed  as 


Ppe  0.8  qO.4 

esuB  =  0-0694  (^)  (4.5) 

^th 

It  is  due  to  the  nozzle  configuration  for  a  submerged  nozzle 
where  S  is  the  submergence  length  in  inches  (the  length  that 
the  nozzle  is  imbedded  in  the  combustion  chamber) . 

6 .  Nozzle  Erosion  Loss  ( ^EROS ) 

This  loss  results  from  nozzle  erosion  during  the 
burn  and  is  expressed  as 


^EROS 


[1 


X  100 

^sPth 


(4.6) 


Ian  is  the  theoretical  ^or  the  initial  nozzle  expan- 

spth  sp 

sion  ratio  and  ^sp^,  is  the  theoretical  Igp  for  the  mean 

expansion  ratio.  ^sp^,  is  determined  by  interpolation  of 

Ian  versus  the  nozzle  expansion  ratio  data  from 

Pth 

"MICROPEP  .  " 


B.  COMBUSTION  EFFICIENCY 

NASA  Report  SP-8064  (Ref.  2)  states  that  the  combustion 
efficienty  (  nc* )  is  determined  by  "-.he  completeness  of  metal 
combustion  within  the  motor  and  by  the  degree  to  which  com¬ 
bustion  products  reach  chemical  eguilibrium  among 


I 


themselves 


r 

» 


t| 


These  effects  depend  primarily  on 


the 


residence  time  (tj-)  as  expressed  in  eguation  2.9.  An  np*  of 


90-95%  can  be  attained  if  residence  time  is  10-15  msec  (see 
Figure  3 ) . 


C.  PROGRAM 

The  subroutines  which  solve  for  the  thrust  and  combus¬ 
tion  efficiencies  are  listed  in  Appendix  A,  Subroutine  PFE 
and  Subroutine  LOSCF.  Basically,  they  followed  the  method 
discassed  in  Chapter  II.  The  algorithm  is  shown  below 

1.  Initialize  Cf  (coefficient  of  thrust)  as  ^f^h 
(theoretical  coefficient  of  thrust  from  "MICROPEP" 
for  shifting  equilibrium) . 

2.  Calculate  A^h  ^th* 


^th  " 


2( 


Cfth^c^ 


1  /2 

) 


3  . 


4  . 


5  . 


6  . 


Calculate  ric^  using  the  relation 

=  1“  (  £div+=bl‘'’^tp+^kin‘*'^sub‘*’^eros 

Calculate  Cf  using 

C  f  —  C  f  Hq 
th  ^f 

Correct  Df^  using 

1/2 

Dth  “  2  - ) 

^f  ’  ^C 

Repeat  Steps  3  to  5  until  Aff,  converges. 


The  reason  that  Af^  used  for  the  iteration  variable  was 

that  He  (through  enr,  e^P '  ^SJIR^  depends  upon  Ai.u. 


V.  GRAINS 


"GRAINS"  is  PC  based  code  for  the  design  of  radial  sided 
star  grains.  The  original  version  was  written  by  G.  J. 

Woten  and  R.  McCortnick  of  CSD,  United  Technologies,  November 

1984.  It  was  then  modified  by  LT  G.  Liston  of  AFWAL/Port, 

December  1984.  The  general  shape  of  this  grain  and  the 

geometric  variables  are  shown  in  Figure  5. 

"GRAINS"  can  be  used  in  either  of  two  modes:  the 

"Design  Mode"  or  the  "Burnback  Model." 

A.  "DESIGN  MODE" 

In  the  "Design  mode"  the  program  solves  for  geometries 
that  satisfy  input  reguirements  on  the  number  of  star  points 
reguired  port  area,  perimeter  factor,  web,  and  etc., 
together  with  the  allowed  tolerances  on  these  parameters. 

The  geometric  constraints  for  the  grain  are  determined 
from  the  propellant  characteristic  and  motor  constraints. 
For  example: 

Abort  =  [|]d2[(i  -  [- - 

^  ■‘•SpFp 

D,  t^  and  F  are  rnotor  inputs  reguirements  and  1,,^  and  Pp  are 
determined  once  the  propellant  has  been  selected. 
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BURNBACK  MODE 


The  "burnback  mode"  uses  the  length  of  the  grain, 
initial  port  area,  perimeter,  and  the  results  from  "Design 
Mode"  to  calculate  the  perimeter,  perimeter  factor,  port 
area,  and  burning  surface  area  as  functions  of  the  web 
burned . 

C.  COMMENTS  ON  PROGRAM  OPERATION 

Care  must  be  used  for  the  inputs  and  input  tolerances, 
otherwise  no  acceptable  solutions  will  be  obtained.  It  is 
best  to  start  with  large  allowable  tolerances  on  Aport'  ©tc . 
(typically  20%).  After  initial  solutions  are  found  in 
"Design"  for  the  number  of  star  points,  range  of  web 
thicknesses,  etc.,  the  range  of  the  input  variables  can  be 
reduced  along  with  the  allowable  tolerances. 


VI.  ROCKET 


"ROCKET"  is  also  a  PC-based  computer  code.  It  was 
originally  written  by  The  Lockheed  Company  for  a  mainframe 
computer.  Later,  it  was  transferred  for  use  on  the  IBM-PC 
by  J.  P.  Francis  of  the  Naval  Weapons  Center  (August  1985). 
The  function  of  "ROCKET"  is  to  calculate  the  pressure-t ime 
and  thrust-time  profiles  for  the  motor,  including  theoreti¬ 
cal  calculations  for  both  ignition  and  tailoff  transients. 
The  description  of  the  motor  can  include  igniters,  multiple 
propellant  grains,  and  motor  insulation.  In  this  thesis 
work  the  program  was  restricted  to  a  motor  with  a  single 
grain  configuration. 

A.  BASIC  ASSUMPTIONS 

The  assumptions  made  in  "ROCKET"  are: 

1.  The  combustion  products  are  ideal  gases. 

2.  Prooellant  burning  rate  follows  the  expression 

.  n 

r  =  aPc  . 


3  . 

There  are 

no 

effects  of 

mass  addition  or  erosive 

burning  in 

the 

combustion 

chamber . 

4  . 

C*  varies 

with 

P(,  in  the 

form  of 

C*  = 

^ref 
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5.  Tae  temperature  sensitivity  used  to  correct 

for  the  effects  of  propellant  temperature  on  the 
propellant  burning  rate  and  chamber  pressure. 

"  ^(70°)  eO  *01  crp^Tp-TO®  ) 

P 

where 

ap  =  (l-n)TrK 

B.  ANALYTICAL  BACKGROUND 

The  program  is  based  upon  the  method  discussed  in 
Chapter  II,  i.e.,  it  numerically  integrates  equation  {2.21) 
to  obtain  ( t )  .  Then  it  calculates  Cf  with  a  fixed  ratio 
of  specific  heats  (^  .  Thrust  is  then  calculated  using 
equation 

F  =  Cf ( t)P^( t)Ath( t) XCp 

where  X  is  nozzle  divergence  loss  and  Cj)  is  the  nozzle 
discharge  coefficient.  X  and  Cp  are  only  two  parts  of 
In  order  to  use  ROCKET  in  the  combined  program,  X  was  set 
to  unity  and  Cp  was  equated  to  rig  • 

C.  INPUTS 

Inputs  to  "ROCKET"  include: 

1.  General  motor  parameters: 

a.  Throat  area  (Af^l-j) 

b.  Nozzle  exit  area  (Ag) 

c.  Total  motor  volume  (Vg) 

d.  Ambient  temperature  and  pressure 


e.  Nozzle  throat  plug  closure  blowout  pressure 

(nozzle  is  oluqged  until  this  pressure  is 

exceeded,  typically  35  psia) 

f.  Initial  pressure  in  the  motor  (Pzero'  typically 
15  psia) 

g.  Nozzle  divergence  loss  X  (set  equal  to  1  in  this 
work ) 

h.  Ratio  of  specific  heats  for  chamber  gases  (y) 

( this  was  obtained  from  "MICROPEP"  chamber 

conditions ) 

i.  Design  pressure  (P^) 

j.  Throat  radial  erosion  rate  (erosr). 

Motor  element  descriptions: 

a.  Propellant  burning  rate  (r) 

b.  Burning  rate  exponent  (n) 

c.  Propellant  temperature  sensitive  (1:^) 

d.  Burning  rate  reference  pressure 

e.  Propellant  C* 

f.  Pressure  correction  exponent  for  C* 

g.  Propellant  density  (pp) 

h.  A  table  of  web  versus  burning  surface  area. 

OUTPUTS 

The  outputs  include: 

A  pressure  versus  time  table 


A  thrust  versus  time  table 


3.  The  free  volume,  propellant  consumed 


( weight ) 


i 

i 


throat  area  versus  time  table 


v' 


V* 


VII.  FLYIT 


"FLYIT"  is  a  three-degree-of-f reedom  (range,  altitude 
and  angle  of  attack)  trajectory  simulation  program  origi¬ 
nally  developed  by  G.  Burgner  of  the  Propulsion  Analysis 
Branch,  Naval  VJeapons  Center  (NWC).  A  User's  guide  was 
later  published  by  Y.G.  Coenen  (Ref.  8).  The  program  can 
simulate  vertical-plane  trajectories  which  include  take-off, 
climbs,  cruises,  landing,  and  rocket  boost.  It  is  a  tool  to 
synthesize  and  analyze  the  trajectories  of  ai r-to-surf ace , 
and  surface-to-surface  missiles. 

The  basic  structure  of  "FLYIT"  consists  of  five  func¬ 
tional  sub-groups,  they  are: 

Atmospherics 
Trajectory  controls 
Aerodynamics 
Rocket  boosters 
Sustanier  propulsion. 

Users  are  allowed  to  select  a  different  option  in  each 
functional  group  to  simulate  different  flight  profiles  with 
different  propulsion  systems.  Through  the  use  of  "FLYIT"  it 
is  possible  to  estimate  the  terminal  velocity  and  range  of 
certain  propulsion  systems. 

"FLYIT"  was  not  incorporated  into  the  main  body  of  the 
present  program  due  to  time  limitations. 
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VIII.  COMBINED  PROGRAM 


The  combined  program  ("SPRMD”)  is  listed  in  Appendix 
and  basically  follows  the  design  process  discussed  in 
Chapter  II.  The  basic  structure  of  this  program  is  shown  on 
Figure  6. 

This  program  combines  "MICROPEP"  (to  calculate  the 
theoretical  thermodynamic  performance  of  the  propellant), 
"PFE,"  and  "LOSCF"  (to  calculate  nc^/  to  estimate 

and  size  the  throat  and  exit  areas),  "GRAINS"  (to  provide 
web  burned  versus  burning  surface  area  for  radially  sided 
star  grains),  and  "ROCKET"  (to  compute  the  pressure-t ime 
and  thrust-time  profiles).  Inputs  to  the  program  include 
"INPUTO"  (the  initial  mission  requirement  inputs),  "INPUTl" 
(inputs  for  "MICROPEP"),  "INPUT2"  (inputs  for  "GRAINS"),  and 
"RCKTIN"  (inputs  for  "ROCKET").  All  of  these  programs  com¬ 
prise  an  interactive  program  for  solid  propellant  rocket 
motor  preliminary  design. 


A.  ALGORITHM 

A  simplified  algorithm  of  this  program  can  be  briefly 
listed  as  follows: 

1.  "INPUTO" 

Input  the  basic  ballistic  and  mission/vehicle 


requirements . 


Calculate  the  required  average 


TERMINAL  VELOCITY 
RANGE 


FIGURE  6. 


FLOW  CHART  FOR  THE  COMBINED  PROGRAM 


requirements.  Calculate  the  required  average 
chamber  pressure  for  a  specific  propellant  aiid 
propellant  temperature. 

"INPUTl" 

Input  propellant  composition,  chamber  pressure,  and 
ambient  pressure.  Create  a  file  for  'MICROPEP." 
"MICROPEP" 

Perform  the  theoretical  performance  calculations  for 
the  selected  propellant. 


Performs  the  following  functions: 

a.  Estimates  the  combustion  efficiency 

b.  Calculates  the  delivered  specific  impulse 

c.  Sizes  the  throat  area,  exist  area,  calculates 
perimeter  factor,  perimeter,  required  web, 
required  burning  rate,  and  etc. 

"INPUT2" 

Input  the  geometrical  description  of  the  star  grain 
and  propellant  characteristics.  Create  a  file  to 
accommodate  the  processing  of  "GRAINS."  Other  grain 
design  programs  (such  as  SPP)  can  be  readily 
incorporated . 

"GRAINS" 


Designs  the  star  grain  and  caJculates  Vi^eb  versus 
burning  surface  area. 


RCKTIN 


Uses  the  results  from  the  previous  programs  to  create 
a  file  for  "ROCKET." 

"ROCKET" 

Calculates  the  pressure-time  and  thrust-time  profiles. 


B.  BRANCHING 

The  algorithm  discussed  above  also  provides  for  the 
capability  of  branching.  For  example,  if  the  desired 
versus  time  is  not  attained,  the  user  can  either  re-do  the 
grain  design  and/or  select  other  propellant  properties. 


IX.  EXAMPLE  OF  DESIGN  PROGRAM 


The  initial  requirements  which  are  input  into  the  pro¬ 
gram  are  shown  in  Table  1,  together  with  a  selected  value 
for  the  temperature  sensitivity  of  the  propellant  (ttk)*  ^ 
reproducibility  tolerance  of  7%  was  chosen. 

The  nominal  operating  pressure  at  Tp  =  70°F  was  calcu¬ 
lated  to  be  1030  psia  and  the  required  propellant  burning 
rate  for  a  web  fraction  of  0.5  was  0.615  in/sec. 

A  composite  propellant  consisting  of  78%  AP ,  12%  HTPB, 
and  10%  Al  was  selected.  "MICROPEP"  was  then  run  for  pres¬ 
sures  of  250  psia  and  the  nominal  P^  value  of  1030  psi. 
Outputs  from  "MICROPEP"  are  presented  in  Table  2. 


TABLE  1 


INPUT  PARAMETERS 


AVERAGE  THRUST,  LBF . 

BURN  TIME,  SEC . 

MINIMUM  THRUST  REQUIRED,  LBF .  8000 

MOTOR  LENGTH,  INCH . 

MOTOR  DIAMETER,  INCH'. . 

MAXIMUM  EXPECTED  OPERATING  PRESSURE... 

TPMAX,  F . 

TPNOM,  F . 

TPMIN,  F . 

TEMPERATURE  SENSITIVITY,  %/F . 

NOMINAL  CHAMBER  PRESSURE,  PSI . 

REQUIRED  BURNING  RATE,  IN/ SEC . 

ALTITUDE,  FT . 

AMBIENT  PRESSURE,  PSI . 

AMBIENT  TEMPERATURE,  R . 


14250.0 

6.50 

.000000 

40.00 

16.00 

1300.0 

150.0 

70.0 

-65.0 

0.200 

1030.2 

0.615 

10000. 

10.13 

483.1 


TABLE  2 


MICROPEP  OUTPUT 
P^  =  10.0  PSIA 


FOR  Pc  =  250  AND  1030  PSIA: 
78%  AP,  12%  HTPB,  10%  Al 


DESIGN 

1  Case 

1  of  2 

OH 

DENS  COMPOSITION 

AMMONIUM 

PERCHLORATE 

(AP) 

-602 

0.07040 

ICL  4H 

IN 

HTPB  (SINCLAIR) 

13 

0.03320 

103H  73C 

10 

ALUMINUM 

(PURE  CRYSTALINE) 

0 

0.09760 

lAL 

INGREDIENT  WEIGHTS  (IN  ORDER)  AND  TOTAL  WEIGHT  (LAST  ITEM  IN  LIST) 

78.0000  12.0000  10.0000  100.0000  , 

INGREDIENT  VOLUME  RATIOS 

AMMONIUM  PERCHLORATE  (AP)  70.487»  HTPB  (SINCLAIR)  22.995» 

ALUMINUM  (PURE  CRYSTALINE)  6.518% 

THE  PROPELLENT  DENSITY  IS  0.06362  LB/CU-IN  OR  1.7610  GM/CC 


NUMBER  OF  GRAM  ATOMS  OF  EACH  ELEMENT  PRESENT  IN  INGREDIENTS 

3.895570  H  0.878965  C  0.663847  N  2.667428  0 

0.370645  AL  0.663847  CL 
CP  «  43.62668610 

aiAMBER  RESULTS  FOLLOW; 

T(K)  T(F)  P(ATM)  P(PSI)  ENTHALPY  ENTROPY  CP/CV  GAS  RT/V  MOL  WT 
3248.  5388.  17.01  250.00  -46.80  240.56  1.1990  3.644  4.667  27.44 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND-  3211.743  AND  3563.987  FT/SEC 

HEAT  CONTENT  (298  REF)  1253.397  CAL/GR  2257.631  BTU/LBM 


0.98321 

H20 

0.71991 

CO 

0.58380 

H2 

0.57339 

HCl 

0.32835 

N2 

0.17947 

A1203* 

0.15900 

C02 

0.11256 

H 

0.07546 

HO 

0.07476 

Cl 

0.00880 

0 

0.00709 

NO 

5.93E-03 

02 

4.51E-03 

AlOCl 

3.65E-03 

A1C12 

3.20E-03 

AlCl 

1.32E-04 

A1C13 

1.28E-04 

C12 

1.19E-04 

AlO 

6.37E-05 

Al 

2.31E-05 

CHO 

1.73E-05 

NH 

1.61E-05 

N 

1.21E-05 

H02 

7.13E-06 

NH3 

5.55E-06 

NHO- 

5.24E-06 

A120 

4.85E-06 

AlH 

2.03E-06 

CNH 

1.91E-06 

A1H02 

1.09E-06 

CH20 

CP  - 

41.92294693 

EXHAUST  RESULTS  FOLLOW: 

i 

T(K)  T(F)  P(ATM) 

P(PSI) 

ENTHALPY 

ENTROPY 

CP/CV 

GAS  RT/V 

2313.  3704.  0.69 

10.13 

-108.91 

240.56 

1.1995 

3.510  0, 

.196 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND-  2650.051  AND  2952.060  FT/SEC 
HEAT  CONTENT  (298  REF)  851.166  CAL/GR  1533.130  BTU/LBM 


I 

i 


1.01958  H20 

0.67108  CO 

0.64916  HCl 

0.59218  H2 

0.33181  N2 

0.20785  C02 

0.18284  A1203& 

0.01785  H 

0.01421  Cl 

0.00497  HO 

0.00233  A1203* 

0.00021  NO 

1.33E-04  A1C12 

1.04E-04  0 

l.OlE-04  AlOCl 

7.51E-05  02 

4.03E-05  AlCl 

1.25E-05  A1C13 

9.06E-06  C12 
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TABLE  2  CONTINUED) 


MICROPEP  OUTPUT  FOR  Pc  =  250  AND  1030  PSIA: 
P^  =  10.0  PSIA.  78%  AP,  12%  HTPB,  10%  A1 


PERFORMANCE:  FROZEN 

ON  FIRST  LINI , 

SHIFTING  ON 

SECOND 

LINE 

ISP 

GAMMA  T* 

P*  C* 

CF 

OPT  EX 

D-ISP 

A*M. 

EX  T 

226.7 

1.1835  2975. 

9.65  5060.1 

}  .44'! 

4.31 

399.2 

0.62924 

1976 

232.5 

1.1497  3045. 

9.77  5082.8 

1.4/2 

4.72 

409.4 

0.63205 

2313 

BOOST  VELOCITIES  FOR 

PROPELLANT  DENSrxY  OF 

0.06362 

(S.G.  OF 

1.T61) 

5./23451. 
60./  7788. 
175. /3i 17. 


10./ia585. 
69./  7128. 
200./3277. 


15./15a57. 
71./  6998. 
300./2336. 


25./12612.  30./11520.  55./  8215. 

88./  6064.  100./  5548.  150./  4113. 
1000./  780.  3000./  269.  5000./  163. 


EXP. 

EXIT 

EXIT 

EXIT 

OPTIMUM 

OPTIMUM 

VACUUM 

VACUUM 

SEA  LV 

SEA  LV 

lATIO 

PRESS 

PRESS 

TEMP 

IMPULSE 

IMPULSE 

IMPULS 

IMPULS 

IMPULS 

IMPULS 

ATM 

SI 

X 

SEC 

SI 

SEC 

SI 

SEC 

SI 

1. 

9.769 

980.6 

3044. 

104.9 

1029. 

195.7 

1919. 

186.4 

1846. 

2. 

3.256 

329.9 

2678. 

175.0 

1716. 

235.5 

2309. 

216.9 

2149. 

3. 

1.215 

123.0 

2386. 

214.4 

2103. 

248.3 

2435. 

220.4 

2184. 

4. 

0.833 

84.4 

2314. 

226.6 

2223. 

257.6 

2526. 

220.5 

2184. 

5. 

0.644 

65.3 

2313. 

234.5 

2300. 

264.5 

2593. 

218.0 

2160. 

6. 

0.523 

53.0 

2311. 

240.7 

2361. 

269.9 

2647. 

214.1 

2121. 

7. 

0.439 

44.5 

2310. 

245.8 

2410. 

274.4 

2690. 

209.3 

2074. 

8. 

0.378 

38.3 

2310. 

250.1 

2453. 

278.2 

2728. 

203.9 

2019. 

9. 

0.331 

33.5 

2309. 

253.8 

2489. 

281.5 

2760. 

197.9 

1960. 

10. 

0.294 

29.8 

2308. 

257.1 

2521. 

284.4 

2789. 

191.5 

1897. 

11. 

0.264 

26.8 

2307. 

260.0 

2550. 

287.0 

2814. 

184.8 

1831. 

12. 

0.240 

24.3 

2307. 

262.6 

2575. 

289.3 

2837. 

177.9 

1762. 

13. 

0.219 

22.2 

2306. 

265.0 

2599. 

291.5 

2858. 

170.7 

1691. 

14. 

0.202 

20.5 

2306. 

267.2 

2620. 

293.5 

2878. 

163.4 

1619. 

15. 

0.137 

19.0 

2305. 

269.2 

2640. 

295.3 

2896. 

155.9 

1545. 

16. 

0.174 

17.7 

2305. 

271.1 

2658. 

297.0 

2912. 

148.3 

1469. 

17. 

0.163 

16.5 

2305. 

272.8 

2675. 

298.5 

2928. 

140.6 

1393. 

18. 

0.153 

15.5 

2304. 

274.4 

2691. 

300.0 

2942. 

132.8 

1316. 

19. 

0.144 

14.6 

2304. 

276.0 

2706. 

301.4 

2956. 

124.9 

1237. 

20. 

0,136 

13.8 

2304. 

277.4 

2720. 

302.7 

2969. 

116.9 

1158. 

21. 

0.129 

13.1 

2303. 

278.8 

2734. 

304 . 0 

2981. 

108.9 

1079. 

22. 

0.123 

12.4 

2303. 

280.1 

2746. 

305.1 

2992. 

100.8 

998. 

23. 

0.117 

11.8 

2303. 

281.3 

2758. 

306.3 

3003. 

92.6 

917. 

24. 

0.112 

11.3 

2302. 

282.5 

2770. 

307;3 

3014. 

84.4 

836. 

25. 

0.107 

10.8 

2302. 

283.6 

2781. 

308.3 

3024. 

76.1 

754. 
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TABLE  2  fCONTINUED) 


MICROPEP  OUTPUT  FOR  P^  =  250  AND  1030  PSIA: 
p  =  10.0  PSIA.  78%  AP,  12%  HTPB,  10%  Al 


DESIGN  1  Casa  2  of  2 


DH 

DENS  COMPOSITION 

AMMONIUM  PERCHLORATE  (AP) 

-602 

0.07040 

ICL  4H 

IN 

40 

HTPB  (SINCLAIR) 

13 

0.03320 

103H  73C 

10 

ALUMINUM  (PURE  CRYSTALINE) 

0 

0.09760 

lAL 

INGREDIENT  WEIGHTS  (IN  ORDER)  AND  TOTAL  WEIGHT  (LAST  ITEM  IN  LIST) 

78.0000  12.0000  10.0000  ,100.0000 

INGREDIENT  VOLUME  RATIOS  ‘ 

AMMONIUM  PERCHLORATE  (AP)  ..  70.487%  HTPB  (SINCLAIR)  .  22.995% 

ALUMINUM  (PURE  CRVSTALINE)  6.518% 

THE  PROPELLENT  DENSITY  IS  0.06362  LB/CU-IN  OR  1.7610  GM/CC 

HUMBER  OF  GRAM  ATOMS  OP  EACH  ELEMENT  PRESENT  IN  INGREDIENTS 

3.895570  H  0.878965  C  0.663847  N  2.667428  O 

0.370645  AL  0.663847  CL 
CP  -  43.75085831 

CHAMBER  RESULTS  FOLLOW: 

T(K)  T(F)  P(ATM)  P(PSI)  ENTHALPY  ENTROPY  CP/CV  GAS  RT/V  MOL  WT 
3375.  5616.  70.08  1030.24  -46.80  230.37  1.1952  3.597  19.485  27,80 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND-  3247.693  AND  3603.495  FT/SEC 

HEAT  CONTENT  (298  REF)  1309.442  CAL/GR  2358.580  BTU/LBM 


1.01084  H20  0.71716  CO  0.59375  HCl  0.57367  H2 

0,32895  N2  0.17941  A1203*  0.16171  C02  0.07526  H 

0.05727  HO  0.05233  Cl  0.00586  NO  0.00506  A1C12 

4.49E-03  0  3.70E-03  AlOCl  3.15E-03  02  2.63E-03  AlCl 

3.14E-04  A1C13  l,83E-04  Cl2  7.24E-05  AlO  5.77E-05  CHO 

3.63E-05  Al  2.75E-05  NH  2.71E-05  NH3  1.55E-05  N 

1.33E-05  H02  9.39E-06  NHO  8.06E-06  CNH  5.08E-06  AlH 

4.45E-06  CH20  3.50E-06  A120  2.16E-06  A1H02 

CP  -  40.18283081 

EXHAUST  RESULTS  FOLLOW: 

« 

T(K)  T(F)  P(ATt'I)  P(PSI)  ENTHALPY  ENTROPY  CP/CV  GAS  RT/V  MOL  WT 
1852.  2874.  0.69  10.13  -130.14  230.37  1.2087  3.491  0.197  28.64 

DAMPED  AND  UNDAMPED  SPEED  OF  SOUND-  2373.392  AND  2644.581  FT/SEC 

HEAT  CONTENT  (298  REF)  659.260  CAL/GR  1137.465  BTU/LBM 


0.98547  H20  0.66303  HCl  0.63198  CO  0.63022  H2 

0.33192  N2  0.24696  C02  0.18531  A12034  0.00097  H 

7.87E-04  Cl  l,10E-04  HO  2.78E-06  A1C12  2.25E-06  NO 

2.21E-06  A1C13  l,52E-06  NH3 


PERFORMANCE:  FROZEN  ON  FIRST  LINE,  SHIFTING  ON 


SECOND  LINE 


TABLE  2  (CONTINUED) 

MICROPEP  OUTPUT  FOR  P^  =  250  AND  1030  PSIA 
=  10.0  PSIA.  78%  AP,  12%  HTPB,  10%  Al 


ISP  GAMMA  T* 
262.5  1.1907  3081. 

269.3  1.1617  3139. 


P*  C* 

39.69  5084.9 
40.09  5126.9 


CF  OPT  EX 
1.661  12.28 
1.690  13.25 


D-ISP  A*M.  EX  T 

462.2  0.15344  1610. 

474.2  0.15471  1852. 


BOOST  VELOCITIES  FOR  PROPELLANT  DENSITY  OF  0.06362  (S.G.  OP  1.7^1) 

5./27165.  10./21528.  15./18368.  25./14609.  30./13344.  55./  9516. 
60./  9021.  69./  8257.  71./  8106.  88./  7024.  100./  6426.  150./  4764. 
175./4224.  200./3796.  300./2705.  1000./  904.  3000./  312.  5000./  188. 


EXP. 

EXIT 

EXIT 

EXIT 

OPTIMUM 

OPTIMUM 

VACUUM 

VACUUM 

SEA  LV 

SEA  LV 

ATIO 

FREDS 

PRESS 

TEMP 

IMPULSE 

IMPULSE 

IMPULS 

IMPULS 

IMPULS 

IMPULS 

ATM 

SI 

K 

SEC 

SI 

SEC 

SI 

SEC 

SI 

1. 

40.059 

4057.9 

3139. 

106.6 

1045. 

197.7 

1939. 

195.4 

1936. 

2. 

13.353 

1352.6 

2722. 

176.8 

1734. 

237.6 

2330. 

233.0 

2308. 

3. 

4.904 

496.8 

2391. 

216.7 

2125. 

250.1 

2453. 

243.3 

2410. 

4. 

3.289 

333.2 

2270. 

229.5 

2250. 

259.4 

2544. 

250.3 

2479. 

5. 

2.435 

246.7 

2183. 

238.3 

2337. 

266.0 

2609. 

254.7 

2523. 

6. 

1.913 

193.8 

2115. 

245.0 

2402. 

271.1 

2658. 

257.4 

2550. 

7. 

1.563 

158.3 

2060. 

250.3 

2454. 

275.1 

2698. 

259.2 

2563. 

8. 

1.314 

133.1 

2014. 

254.6 

2497. 

278.5 

2731. 

260.3 

2579. 

9. 

1.129 

114.3 

1975. 

258.3 

2532. 

281.4 

2759. 

260.9 

2584. 

10. 

0.986 

99.9 

1940. 

261.4 

2563. 

283.8 

2783. 

261.1 

2586. 

11. 

0.873 

88.4 

1910. 

264.2 

2590. 

286.0 

2805. 

261.0 

2585. 

12. 

0.781 

79.1 

1882. 

266.6 

2615. 

287.9 

2824. 

260.7 

2582. 

13. 

0.706 

71.5 

1858. 

268.8 

2636. 

289.7 

2341. 

260.1 

2577. 

14. 

0.643 

65.1 

1835. 

270.8 

2655. 

291.3 

2856. 

259.4 

2570. 

15. 

0.589 

59.7 

1814. 

272.6 

2673. 

292.7 

2870. 

258.6 

2562. 

16. 

0.543 

55.0 

1795. 

274.3 

2690. 

294.0 

2883. 

257.6 

2552. 

17. 

0.503 

51.0 

1778. 

275.8 

2705. 

295.2 

2895. 

256.6 

2542. 

18. 

0.468 

47.4 

1761. 

277.2 

2718. 

296.4 

2906. 

255.5 

2531, 

19. 

0.438 

44.3 

1746. 

278.5 

2731. 

297.4 

2917. 

254.2 

2519. 

20. 

0.411 

41.6 

1731. 

279.8 

2743. 

298.4 

2927. 

253.0 

2506. 

21. 

0.386 

39.1 

1717. 

280.9 

2755. 

299.4 

2936. 

251.6 

2493. 

22. 

0.365 

36.9 

1704. 

282.0 

2766. 

300.3 

2944. 

250.2 

2479. 

23. 

0.345 

34.9 

1692. 

283.0 

2776. 

301.1 

'  2952. 

248.8 

2465. 

24. 

0.327 

■33.1 

1681. 

284.0 

2785. 

301.9 

2960. 

247,3 

2450. 

25. 

0.311 

31.5 

1670. 

284.9 

2794. 

302.6 

2968. 

245.8 

2435. 

•' 


JOT?? 


The  output  from  "MICROPEP  indicated  the  following 

C*  =  5126.9  ft/sec 
th 

=  269.3  Ibf  sec/lbm 

®Pth 

e  =  13.25 

Moles  of  condensed  AI2O3/IOO  gni  =  0.1794 
Pp  =  .0636  Ibm/in^ 

The  AGARD-AR-230  empirical  equations  were  then  used  to 
estimate  the  value  for  nq  *  chosen  to  be  0.93.  For 

c  the  calculated  residence  time  of  7.2  msec,  Figure  3,  indi¬ 
cates  that  a  better  value  may  have  been  approximately  0.87. 
As  discussed  above,  the  losses  are  iterated  with  the  throat 
area.  The  resulting  losses  and  throat  area  are  presented  in 
Table  3. 

Once  the  losses  have  been  estimated  Igp  can  be  deter¬ 
mined  from 

Isp  =  Isp^^ 

The  following  parameters  were  then  calculated: 

™  =  F/Igp 

V  z  =  Vp/V3 
Dpe  = 

J  =  Ath/Ap 

These  parameters  are  shown  in  Table  4. 
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TABLE  3 


ESTIMATED  PERFORMANCE  LOSSES 


PERCENT  DIVERGENCE  LOSS . 

PERCENT  TWO  PHASE  FLOW  LOSS . 

PERCENT  BOUNDARY  LAYER  LOSS . 

PERCENT  KINETIC  LOSS . 

PERCENT  SUBMERGENCE  LOSS . 

PERCENT  EROSION  LOSS . 

THRUST  COEFFICIENT  EFFICIENCY . 

Etacf=  0.924  Cfc=  1.562  Ath= 

PERCENT  DIVERGENCE  LOSS . 

PERCENT  TWO  PHASE  FLOW  LOSS . 

PERCENT  BOUNDARY  LAYER  LOSS . 

PERCENT  KINETIC  LOSS . . 

PERCENT  SUBMERGENCE  LOSS . 

PERCENT  EROSION  LOSS . 

THRUST  COEFFICIENT  EFFICIENCY . 


1.704 

3.202 

1.828 

0.849 

0.000 

0.007 

0.924 

8.856 

1.704 

3.199 

1.783 

0.849 

0.000 

0.007 

0.925 
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TABLE  4 


CALCULATED  MOTOR  DATA 


THRUST  COEFFICIENT . . 

ETACF . 

SPECIFIC  IMPULSE,  LBF  SEC/LBM . 

MASS  FLOW  RATE,  LBM/SEC . 

PROPELLANT  VOLUME,  CU  IN . 

NOZZLE  THROAT  AREA,  SQ  IN . 

NOZZLE  EXIT  AREA,  SQ  IN . 

EQUIVALENT  PORT  DIAMETER,  IN . 

PORT  AREA,  SQ  IN  . 

PERIMETER  FACTOR . 

REQUIRED  WEB,  IN . 

REQUIRED  WEB  FRACTION . 

VOLUMETRIC  LOADING . 

LENGTH  TO  DIAMETER  RATIO . 

RESIDENCE  TIME,  MSEC . 

THROAT  TO  PORT  AREA  RATIO  (J  FACTOR) . 
REQUIRED  BURNING  RATE,  IN/SEC . 


1.563 
0.925 
231.587 
61.532 
6286.8 
8.851 
117.290 
7.476 
43. ‘893 
2.489 
2.500 
0.312 
0.782 
2.500 
7.256 
0.202 
0.385 


>  .*  J'  -•*  -•  ^  ^  ^  -»  -•  -*  ^ 


At  this  point  the  following  expressions  remain  to  be 
satisfied : 

Perimeter  Factor  =  PFAC  =  Pwetted^ ’^^pe 

‘PFAC  =Ai3/TTDpgL 

"  ^wetted  ^  ■ 

™  “  Pp  ^b 

^  cnom 

r  =  web/^5 

In  these  three  expressions  the  unspecified  variables  are 
PFAC/  A5,  ir,  and  web. 

For  a  selected  web  or  r  the  other  parameters  can  be 
determined.  For  example,  in  the  current  problem  Table  E 
summarizes  the  required  values  for  selected  web 
thicknesses . 


TABLE 

_5 

REQUIRED  PARAMETERS 

FOR  SELECTED  WEBS 

web ( in ) 

r  (in/sec) 

Ak  (in^) 

PFAC 

3.5 

.538 

1796 

1.91 

3.0 

.461 

2095 

2.22 

2.5 

.385 

2514 

2.67 

Once  the  required  burning  rate  (r)  is  determined,  appro¬ 
priate  values  of  a  and  n  are  required  for  the 

.  n 

expression  r  =  ^P^  . 

The  grain  design  problem  is  to  find  a  geometry  which 
will  satisfy  Table  5,  together  with  the  required  conditions 
on  L,  D,  Pf,-t  and  F-t  (and  stress). 


Ss 


In  the  present  program,  only  radially  spoked  star  grains 
can  be  analyzed.  It  may  not  be  possible  to  satisfy  the 
motor  requirements  with  this  type  of  grain.  The  process  is 
essentially  one  of  tr ial-and-error . 

Three  attempts  at  meeting  the  design  requirements  are 
presented.  One  6-point  and  two  8-point  designs  were  tried. 

A  typical  output  from  the  "Design  Mode"  of  "GRAINS"  is 
shown  in  Table  6. 

Once  the  "GRAINS  Design  Mode"  has  been  run,  "GRAINS 
Burnback"  is  run  to  obtain  A5  versus  web. 

The  output  from  "Burnback"  for  Case  II  is  shown  in  Table 


"ROCKET"  can  then  be  run  to  obtain  the  pressure-t ime  and 
thrust-time  profiles.  Input  data  for  Case  II  are  shown  in 
Table  8. 

A  small  section  of  the  output  from  "ROCKET"  is  shown  in 
Table  9  for  Case  II. 

A  summary  of  the  results  for  the  three  cases  is  pre¬ 
sented  in  Table  10  and  pressure- t ime  traces  are  shown  in 
Figure  7. 

None  of  the  cases  attempted  met  the  desired  design 
values  with  sufficient  accuracy.  Further  iterations  may  re¬ 
sult  in  a  better  solution.  However,  the  star  design  has 
considerable  sliver  with  the  resulting  long  tailoff.  This 
section  was  presented  in  order  to  demonstrate  the  process 
required  to  obtain  a  final  grain  design. 


TABLE  6 


%  N.  •-'  *-  O  O 


OUTPUT  FROM  "DESIGN  MODE"  OF  "GRAINS" 


Uesis.-:  . 

::DUtS 

*aX13UD 

Nuade'  of  5b 

ekes  =  5 

brain  Outer  Radius  =  9. 

Desired 

Port  Area  = 

43.89293 

Desired  Periieter  Factor  = 

RIl  = 

3.S00GC 

DEL2  =  ■ 

-0.25000 

FIN2  =  2.00000 

SRIl  = 

0.25000 

DELI  = 

C. 10000 

FINi  =  0.25000 

SRI2  = 

0.25000 

DEL4  = 

0.10004 

FIN4  =  0  25000 

FINHEB  = 

4.00000 

yEBSTP-= 

0.25000 

KSFQKE  =  0.50000 

CIRTQL  = 

0. 10000 

PTTGl  = 

0. 10000 

“cRIQu  -  v.^OOOO 

Icier  ar.: 

e  .lilts 

“or:  lirea:  r!in:»uii  = 

39.50364 

.laxiiusi 

Feriiete 

r:  Itiniiui  = 

52.61331 

Raxiaui 

=  64.20515 

Circle 

:  ^iniiua  = 

324.00000 

r.axiBuo 

=  395.00000 

2.48915 


•V-  V  ■  V  -- 


TABLE  6  (CONTINUED) 

OUTPUT  FROM  "DESIGN  MODE"  OF  "GRAINS" 


H  Afcr  -ERI!*  "AC  RRl  R2  R3 

ERROR  CC  zrti'.i’  tnar;  FERTG.,  :r!r’'jfcrr  ecual. 
3  -3. AES  4;. 131  2.01o  3.500  3.750  5.00C  5.150 
ERROA  Cl  create'  tfiar.  FERTG...  t.'ie'efc'e  ecaal 


iiel  3R. 


CDk"^  Sh'V  ?U' 


-.BEL  -'.E;- 


ERROR  ->  =ci4 
E3u»:fi:  E  ' 
ERROR  =h!4 
Ecoees:  E  F 
ERROR  -A  F.C2  4 
icoses:  S  P 
ERROR  Pci  4 
ssotea:  3  ? 
Pr'i4 

Soofc^:s:  a  F 


.  0.0.  less  tnar  u.O 

il:  0.!  Phil:  0.3  rhi3:  0.1 

1  O.Ci  ieSs  tnan  0.0 

til  u..  Pcil;  C.i  'Ct3;  0.1 

3  0.01  less  tf'ac  0.0 

il:  0.1  ?Ci2:  0.5  ?ci3:  j.l 

I  -O.li  less  tnan  0.0 

I I  *  0 « *  •  "'  •  *  ^  ri  i  3  j  0 » 1 

1  -C.i;  \iii  tear  O.C 

il:  0.1  Pci2:  0.5  =r:3:  O.i 


TABLE  7 

GRAINS  "BURNBACK"  OUTPUT  FOR  CASE  II 


-t5;;r  . 

ncuis 

haxuun 

Nusoer  o' 

SpC-ncS  -  0 

jrain 

Qutef  Raciu 

5  =  S.C'OCCO 

Sesired 

Port  Area 

=  i:.:::oo 

Des.r 

e;  Periaeier 

.‘actor  =  2.553CC 

RIl  = 

3.00000 

3EL2  = 

0.  N/OOOO 

r  I  i’ll  = 

O.OOGOO 

SSI;  = 

0.25000 

del:  = 

0.00000 

.‘INI  = 

C. 25000 

CS TO  s 

C. 25000 

DEL4  = 

0.00000 

riN4  = 

C.25COO 

"INitEB  = 

0.00000 

kiESSTF  = 

0.00000 

KSPC.<£  = 

O.COOOO 

CIRTCi.  = 

j. 00000 

P-TCL  = 

0.00000 

PERTGl  = 

0.00000 

Surnoack 

InButs 

• 

train  Lenqtn  =  ■ 

tO. 00000 

PHI2  =  2S.S1A0C 

?HI4  =  2. 30300 

STEPl  = 

e 

STEP2  =  5 

QTC 

t  M  i 

<5  -  c 

TTrC-*,  z 

Nuiber  of  Spokes  =  8 

•ei  Thickness  =  2.5C000 

Web  Perneter  Periieter  Burn  Port 
Burned  Factor  Area  Area 

(IN)  (IN)  (SQ  IN)  (SJ 


Burnback  of  SRI  (Spone  Corner  Radiui) 


0.000 

58.737 

2.552 

2349.472 

42.121 

0.050 

53.051 

2.481 

2362.038 

*rB.  ObO 

0.100 

53.365 

2.416 

2274.604 

48.027 

0.15C 

53.679 

"  ■’57 

2387.170 

51. 003 

0.200 

59.393 

2.303 

2399.737 

53.334 

0.250 

60.308 

2.253 

2412.303 

57.002 

0.250 

60.308 

2.153 

2412.303 

57.002 

0.465 

61.383 

;.C54 

2455.576 

71.054 

0.673 

£2.150 

1.904 

2425.997 

34.747 

0.634 

62.589 

1.784 

2503. 564 

37.373 

1.103 

1.682 

2508.279 

1 1  O'.  643 

»  •  wa.w 

61.504 

■ 

2500. 14! 

:;:.£33 

60 


'>•  y  y  y  y 
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TABLE  7  (CONTINUED) 

GRAINS  "BURNBACK"  OUTPUT  FOR  CASE  II 


ifnDacx  c;  ;;:cke  iica.- 

..323  02.504  :.57S  2500.141  124. 85£ 


447  2364.630  13:. SI? 
366  2255.950  136.593 


cafnOdC  k 

of  raiiiainae'  to  «£ 

1  -  C 

I.cEq 

"•4 

Z  •  ^C'l: 

*  •  w  *  4 

1.848 

52.036  1.1B9 

202l.4S'3 

2.011 

A*?*?  ' 

2061. 309 

2.174 

51.566  1.119 

2053.440 

*  •  vw  4 

ff'  '  A'iC 

M44tfWV  «4V«« 

2C77.136 

4  4  i'-JM 

«Jte4^t3  •4VW. 

2;i7. 374 

I*ro 

ThiOtf’rSS  . 

Percen: 

0!  fual  rriaininc  = 

9.40  , 

i 

BurnDaci 

■jf  resaininc  !uf. 

2C.320  :.i05 

0*4,0  i 3o 

*»  Q/*' 

3*»w 

•Z  '-r.  r 

^  0  Wv*  3 

0*7 0 154 

3. 066 

10.419  C,2C'S 

•  '  '  TC  * 
•*W4  /.JJ 

3.200 

5.272  O.IOS 

A  4  .  «  *.  T 

«  <  n 

•*•  A'V.*.  •.  AAA 

/-  A  •• 

v*  V . 


TABLE  8 

INPUT  DATA  FOR  "ROCKET",  CASE  II 
•■otcr  Ptfrfcrsanc?  -rograh  <ersi:-r  1., 

DAin..< 

•-"t:"  "AS  :  SPAINS 
SESCPIFTIuN  Cf  uPAIN  I 

BilRN  RATE  T%3346 

BURN  RATE  EKP.  0.40CC 

PI  SUB  K  C-2000 

BURN  RE*.  PRESS.  1030.2 

;  ;TAR  4763.0 

C  STA®  E7?.  3.C061 

3EN3:T\  0.0E33 

ISNI'ICN  TI*E  0.0000 


4763.0 

0.C061 

0.G633 

C.OOOO 

V  .  VvOO 


"E:.TA  IcN.  TII^E  v.vv'OU 

Noraai  Incut  la  f:nisned 
Seneral  Ccnficuratio.n  PARA’IETEP.S  are  ... 


'n'oa;  AREA 

3.3510 

Esr.  AREA 

117.2903 

Eicansion  RATIO 

•  nc*^ 

Total  flotor  VOLUME 

3042.5 

AMSIES*  Te»aera:ure 

70.  V 

AMEIEN*  "reiSure 

* 

Closure  t.udQUT 

35.0 

®:ero 

15. C 

LaBCga 

1.0000 

Cd 

C.5246 

lasisa 

*  «  *  " 

Tnroat  DESISN  -'essur 

.  A  A  ^ 

’nr(0AT  cROS.  HE. A/  TI 

iOOv. 000 

RSS3URE  VS.  THROA"  EROSION 

RATE 

600. 

0.000000 

i400. 

O.CCCOOC 

yEB 

BURN  AREA 

0.000 

2349.472 

0.250 

2411.952 

1.108 

2507.744 

1.540 

2296.212 

2. 011 

2060.360 

2,500 

2037.107 

2.832 

812.473 

3.342 

0,000 

"OTA;.  PROPELLANT  WEI6HT= 


4:1.553  I6NITER  WE16HT. 


.  »  •  •  s.”  o  •  •  -  •  •  .  "  .  *  "  -  -  ^  •  »  •  *  •  .  -  •  "  v.*  -  '  •  V  •  *  -v"  •  «  "  •  •  ‘  '  -v 


TABLE  9 


SAMPLE 

OUTPUT 

-FROM  " 

ROCKET" 

,  CASE 

II 

Tifr-: 

r* 

PCC7 

r i  i  V  u  L 

-  iiEIurT 

hIZj"-*' 

if'w  ft 

ilVJ  . 

A 

AVfl.. 

7.;C5 

536.14 

572.85 

-655. '.3 

-  « c 

379.60 

33.81 

0.00 

OP  n 

3.851 

0.0000 

7.225 

-545.45 

7543.30 

nq.A 

33.04 

0.00 

j4.  JO 

8.351 

0.0000 

7.246 

553.73 

-634.70 

7554.42 

380.95 

32.27 

0.00 

At  er 

JO,  JJ 

8.851 

0.0000 

7.266 

546.35 

-624.11 

7564.68 

381.60 

31.52 

o.co 

22.78 

8.BS1 

0.0000 

7.287 

534.18 

-613.15 

7574.71 

332.24 

30.78 

0.00 

32.03 

8.851 

0.0000 

7.3C7 

521.73 

-602.31  . 

:’584.50 

282.56 

3C.  06 

(,,00 

A  ♦  AQ 

0,  • 

8.851 

0.0000 

7.327 

503.51 

-5^1  2^ 

7534.07 

383.47 

29.35 

0.00 

Jv,  Ob 

a.S51 

c.oooo 

-^iO 

/  ■  JTu 

437.51 

C  'N  ^  *. 

« ■  c ^ 

7603,40 

384.07 

28.55 

0.00 

29.84 

S.351 

0.0000 

7.36S 

.•nc 

'V  W  W  «  /  j 

-569.35 

/  6 « 4!  ■  5  X 

384.55 

27.37 

0.0c 

23.14 

3,85; 

V, JuOv 

7.388 

“74,  17 

-559.22 

^  '  -A  .  it 

''.;c  Af 

27,30 

0.00 

z 

.  3.851 

•  ••.AAA 
V,  VWV 

7.40S 

462.32 

-548.63 

7630.03 

^  w 

26.64 

u.  00 

27.77 

S.351 

0.0000 

7.430 

451.69 

-535.30 

7636.56 

ygSi, 

-AC  uO 

0.00 

27,11 

6.851 

l.COCO 

-•  lC•^ 

.  .  --y 

44;. 77 

-52i. 01 

764o. S3 

wuO , 

Ow 

0.0? 

21.4c 

2.251 

/•  A.'-  • 

V  ,  ,  V  V  V 

7.  4/  i 

“  JVl  ‘v'O 

-317.30 

';54.S5 

00/ ♦ 

A  T  •  - 

0,  viz 

teW  1  Cy 

;•  .■^c  • 

,  8  J  i 

V  . 

7. '31 

4 1  3  *  mO 

*507. S£ 

>'  6  u  y  •  /  6 

,C** 

A  >  «  ‘A 

i*T,  ;,i 

0.00 

25.13 

*'  AC  ' 

8,  CJi 

■  J  V  W 

••  C<  0 

7  •  J.  4 

-•i  55, 0** 

7670.42 

383.33 

•  wO 

C,  vO 

24.53 

2.851 

o.ooco 

-»  C'?'> 

/  • 

233.16 

-453.28. 

7677.31 

333.81 

22.95 

o.-co 

22.97 

3.851 

A  AAA.  A 

V, vvvv 

7.u73 

0/ 

-469.21 

tA 

CtS^tfl 

2l,£l 

C.  00 

22.30 

^  AC* 

u ,  00, 

V, vvVV 

7.cil4 

360.71 

ve.' 

*4jv.  /V 

7706.01 

390.59 

AA  tAl 

4V, 

O.oO 

21.68' 

3.351 

•..0000 

7 «  um^ 

342.62 

■“Z2, 72 

7713.01 

291.42 

19.63 

0.00 

20.60 

3.351 

0.0000 

7.636 

225.25 

-415.25 

7721.25 

13.67 

0.00 

19.56 

S.351 

O.O'OOO 

7.737 

308.59 

•4:5. 

7743.06 

232.35 

♦  t 

././I 

0.00 

12.56 

3 ,  d5 i 

.  0000 

7.778 

232.52 

-3;.. 39 

*•  t  e  .•  <  c 

4O 

293.55 

1 :  **a 

•  w>  /  9 

0,00 

17. d1 

S.35i 

V ,  V’Jvv 

7.813 

277.30 

*365.36 

7764.58 

334.33 

15.90 

A  A,-"' 

V*  VV 

16.63 

6.351 

V,  'oOOO 

7.860 

262.63 

-350.52 

7774.54 

234.55 

15.05 

o.co 

15.81 

8.231 

’.301 

248.58 

-335. 56 

7794.07 

3:5.56 

!"•*.  24 

0.00 

14,':7 

2.851 

■■  AAA,'- 

vVVv 

7.342 

>4*?^  »  0 

*ttlj  «  «  w 

-321.03 

^792.33 

12.47 

0.00 

14, 17 

A  AT  • 

Z% 

V  . 

7.333 

5.024 

''At  AT 

*i<;/  •  VO 

*«  Ty • JV 

'6v'..42 

396.65 

2:7.17 

» ■  tA 

0.00  ' 

*  •  *  /\ 

•  AC< 

0,  CJi 

A  ,*>P  • 

S.COi 

•A  •  A  '• 

203.57 

750;, 35 

::.oi 

o.co 

12, So 

0.0000 

9.065 

195.22 

-280.38 

’316.30 

337. o5 

11.34 

0.00 

11.96 

2. 851 

C.^IOOO 

3.106 

186.33 

-267.70 

7323.99 

•  296.10 

Iv/.  69 

0.00 

•  *  AQ 

«  ,  •  ,>0 

8.B51 

•..■9000 

S.  147 

176.23 

-253.44 

/  2wV • 6  ^ 

W  Ju  •  W 

l'v,C7 

0.00 

10.64 

6.  B5l 

(..OCOO 

S.1S8 

165.21. 

-221.51 

7536.97 

*5C|0  O') 

w  ^ 

9.54 

c.oo 

10.04 

8.351 

C, vCvO 

5.243 

155.54 

-133.32 

7845.36 

395.50 

•’  9.09 

o.co 

3.40 

B.851 

c.oooo 

5.341 

145.01 

-98.55 

7552.72 

400,31 

'  3.55 

0.00 

8.76 

8.851 

0.0000 

8.464 

133.96 

-84.16 

7874.60 

401.32  - 

7.31 

0.00 

b.lO 

8.851 

0.0000 

8.628 

120.86 

-76.13 

7393.38 

4C2.55 

7.14 

0.00 

7.31 

8.851 

c.oooo 

8.732 

103.37 

-63.82 

7911.45 

402.65 

6.43 

0.00 

•  p/* 

8.651 

o.cooo 

8.956 

37.33 

-64.30 

7327.17 

404,66 

5.78 

0.00 

5.33 

3.851 

c.oooo 

; .  11  '3 

87.73 

-59.36 

7941.29 

405.56 

c 

'J»  k  i 

c.oo 

c  aa 
y ,  Jy 

2.351 

0.0000 

3.324 

76.36 

-52.70 

7356.39 

406.55 

•  C  i 

*  •  ^  i 

0,00 

4,  £j 

B.35; 

C .  O^O'Cu 

3.523 

66.16 

-4’. 05 

7970.42 

407.41 

3.31 

•..00 

4,02 

2.  Scil 

0,  CvCO 

3.775 

55.33 

"•»v.  CO 

7964.27 

4  j5. jO 

~  At 

0.00 

9.27 

6.851 

7. 

10.062 

44.65 

•  -24.30 

i;:-  cp 

.  ^  /  I  ww 

40:. 14 

-.54 

■0.00 

-.72 

8.85. 

. .  I'COC 

10.383 

34.48 

-27.31 

8009.53 

4C9.90 

2,.  O-i 

0.00 

^  1 

3,2-:: 

^ , 000c 

10.733 

74. 48 

» t  lo 

6'.'.; , .  65 

4;-.6: 

•  r 

0.00 

•  PO 

A  t  jJ 

2.251 

.  .  vv'-.' 

11. 311 

15.62 

-1..34 

5C3v.OO 

-  1  .  .  1  7 

,  3  1 

0.00 

■i .  3“ 

-  ->c  . 

«  •  w  J  . 

«  .  Vv’V  ! 

63 
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TABLE  9  (CONTINUED) 


SAMPLE  OUTPUT  FROM  "ROCKET" ,  CASE 


P=E3’3-- 

>r;Cr 

RACrO?! 

4A3 

TiiIl 

PC 

■1  “ 

rz 

r 

7.573 

379.55 

«  on 

4368.43 

7.&14 

360.71 

3,30 

41 1 1. 61 

7.655 

.jn^i 

2*24 

3373.83 

7.656 

w>  2  < 

3646.54 

7.737 

303.53 

2*41 

3-1';.  23 

7. 773 

^  •  c«. 

0  ‘C 

2i  *Tm 

3:: 1.64 

".515 

r^’TZT’ 

J «  'S 
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SUMMARY  OF  RESULTS  FOR  THREE  GRAIN  DESIGNS 
_ D  =  16  IN,  L  =  40  IN _ 


GRAIN  3 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

An  interactive  computer  code  for  the  preliminary  design 
of  solid  rocket  propellant  rocket  motors  ("SPRMD")  was  suc¬ 
cessfully  developed  and  its  use  was  demonstrated  through  a 
design  example.  The  code  incorporated  several  existing  PC 
codes  and  AGARD  recommended  empirical  expressions  for  losses 
in  aluminized  propellant  systems. 

B.  RECOMMENDATIONS 

Further  modifications  to  improve  the  capabilities  of  the 
code  are  essential  in  order  to  make  "SPRMD"  a  more  effective 
preliminary  design  tool. 

There  is  only  one  grain  configuration  (radially  spoked 
star  grains)  incorporated  into  "SPRMD." 

It  is  necessary  to  add  other  options  for  grain  design 
in  order  to  increase  the  grain  design  flexibility.  One 
possibility  would  be  to  incorporate  the  Grain  Design 
and  Internal  Ballistics  (GDB)  module  of  Solid  Per¬ 
formance  Prediction  Code  (SPP)  (Ref.  3). 

"SPRMD"  also  does  not  currently  incorporate  any  plot¬ 
ting  capabilities.  There  is  definitely  a  need  for  this 
program  to  have  plotting  capabilities  for  the  grain 
design  results  and  the  thrust  and  pressure  versus 
time  profiles. 


"FLYIT"  was  not  incorporated  into  the  program,  although 
"SPRflD"  was  structured  to  have  flight  simulation  capa¬ 
bilities.  It  is  also  recommended  that  "FLYIT"  be 
incorporated  into  "SPRflD"  as  another  enhancement. 
"SPRMD"  is  a  primative  type  of  preliminary  solid  pro¬ 
pellant  rocket  motor  design  code.  Further  enhance¬ 
ments,  such  as  propellant  searching,  screen  editing, 
and  increased  solution  speed  for  "MICROPRP"  and 
"ROCKET"  would  make  this  code  more  user  friendly. 
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program  SPRMD- 

c*********************************************************** 

c  * 

c  SOLID  PROF3LLANT  ROCKET  MOTOR  DESIGN  * 

C  THIS  PROGRAM  (SFRMD)  SOLVE  THE  PRELIMINARY  DESIGN  PROBLEM* 
c  OF  SOLID  IROPELLANT  ROCKET  MOTOR.  SUBROUTINES  INCLUDES:  * 

C  (1)  INPUTO:  BASIC  INPUT  FOR  DESIGN  * 

c  (2)  INPUTl:  CHEMICAL  COMPOSITION  INPUT  * 

c  (3)  MICROPEP:  PERFORME  THE  PROPELLANT  EVALUTION  * 

C  (4)  PFE:  SIZE  THE  THROAT  OF  NOZZLE  PORT  AREA  OF  GRAIN  ETC* 
c  (5)  LOSCF;  CALCULATE  THE  THRUST  COEFFICIENT  * 

c  (6)  INPUT2;  STAR  GRAIN  INPUT  * 

c  (7)  GRAINS:  PERFORM  THE  SPOKED  STAR  GRAINS  DESIGN  * 

C  (8)  RCKTIN:  input  THE  DATA  FOR  THRUST  AND  PRESSURE  PRO-  *  • 

G  FILE  CALCULATION  * 

c  (9)  ROCKET:  COMPUTED  THE  THRUST  AND  PRESSURE  PROFILE  * 

c*********************************************************** 
c 

real  f avg , tb , 1 , d , meop , pc , rrqd , pamb , tamb , h , Id , webrqd , wf , m 
common  /micrp/  a (12 , 12) ,kr(20) ,amat (10, 12) , jat ( 12) , in, is, 
lfie(10,6) ,ie(10,6) , alp (12) ,w27,n,dh(10) ,rho(lO) , 

2wate(10) ,wl(6) ,w43 , ig,np,vnt(201) ,w47,name,ser, 

3floor, itag(lOO) ,wing(10) 
c  integer 

character*10  yea 

c*********************************************************** 

c  * 

c  BASIC  INPUT  FOR  MOTOR  DESIGN  * 

c  * 

C********************************************* ******** ****** 

1000  call  input 0 ( f avg , tb , 1 , d , meop , pc , rrqd , pamb , tauub , h , f min , tpmax , tpnom , 


*tpmin,pik) 
write(*,l) favg 

1  format (lx, 'AVERAGE  THRUST,  LBF .  ’,fl2.1) 

write ( * , 2 ) tb 

2  format (lx, 'BURN  TIME,  SEC .  ',fl2.2) 

write(*, 12) fmin 

12  format (lx, 'MINIMUM  THRUST  REQUIRED,  LBF .  ',fl2.6) 

write(*, 3) 1 

3  format  (lx, 'MOTOR  LENGTH,  INCH .  ',fl2.2) 

write  (* , 4 )  d 

4  format  (lx, 'MOTOR  DIAMETER,  INCH .  ',fl2.2) 

write (*, 5) meop 

5  format! lx, 'MAXIMUM  EXPECTED  OPERATING  PRESSURE. . .  ',fl2.1) 
write ( * , 13 ) tpmax 

13  format  (lx, 'TPMAX,  F . .  . . ',fl2.1) 

write (*,131) tpnom  * 

131  format (lx, 'TPNOM,  F . • .  ',fl2.1) 

write (*,132) tpmin 

132  format (lx, 'TPMIN,  F .  ',fl2.1) 

write(*,14)pik 

14  format (lx, 'TEMPERATURE  SENSITIVITY,  l/F .  ',fl2.3) 

write (* , 6) pc 

6  format (lx, 'NOMINAL  CHAMBER  PRESSURE,  PSI .  ',fl2.1) 

write (*, 7) rrqd 

7  format (lx, 'REQUIRED  BURNING  RATE,  IN/SEC .  ',fl2.3) 

write(*,8)h 

8  format (lx, 'ALTITUDE,  FT .  ',fl2.0) 

write ( * , 9) pamb 

9  format  (lx, 'AMBIENT  PRESSURE,  PSI . ',fl2.2) 
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write (*, 11) tamb 

11  format  (lx, 'AMBIENT  TEMPERATURE,  R . •,fl2.1) 

write (*,10) 

10  format (lx, 'ANY  CORRECTION  OF  BASIC  INPUT?  Y-YES,  N-NO...  ') 

read (*,20) yes 
20  format(alO) 

if (yes.eq. 'y' .or.yes.eq. 'Y')  go  to  1000 
c 

C******** ******** ******************************************* 

c  * 

c  INPUT  THE  DATA  NEEDED  BY  MICROPEP  * 

C  * 

C* *************************** ******************************* 

C 

c 

1010  call  input 1 (pc, pamb) 
write(*,40) 

40  format (lx, 'Are  there  any  corrections  of  input  for  micropep?,y-yes 
*n-no  ' ) 
read ( * , 20) yes 

if (yes. eq. ' Y' .or .yes. eq. 'y ' )  go  to  1010 
c 

c*********************************************************** 

c  * 

c  COMPUTE  THERMAL  EQUILIBRIUM,  PROPELLANT  CHARACTERISTIC,  * 
c  AND  ETC. (USING  MICROPEP,  UNTILL  GOOD  RESULTS  ARE  REACHED)* 

C  * 

C*********************************************************** 


call  micropep 
write(*,50) 

50  formatdx, 'Are  the  results  good?  y-yes  n-no  ') 
read(*,20)ye3 

if (yes.eq. 'y ' .or.yes.eq. ' Y • )  then 
go  to  1020 
else 

•itef*.60) 

Change  the  Basic  Inputs  of  motor.')' 


55 

write (* , 60) 

60 

format(lx, ' (1) 
write(*,70) 

Change 

70 

format(lx, ' (2) 
write (*,80) 

Change 

80 

format (lx, ' 

type  1 

write (*,90) 
format (lx, ' 
read(*, 100) i 

90 

that  M 

100 

format (il) 

if (i.ne. 1. and. i < 

.  ne.2) 

write(*,*)'  Error  inpu^f!  1,  or  2  must  be  typed' 

go  to  55 
el'a 

if(i.eq.l)  go  to  1000 
if(l.eq.2)  go  to  1010 
endif 
endif 

c*********************************************************** 

c  * 

C  CALCULATE  THE  LOSS  PERFORMANCE  &  COMPUTE  THE  DATA  FOR  * 
c  GRAIN  DESIGN  * 
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1020  call  pfe(favg,tb,l,d,pc,cfc,ath,atacf,m,vp,vl,dpe,aport, 
*pf ao , isp , cf th , webrqd , wf , tc , rrqd , ae , atac , xx) 
write (*,50) 
read(*,20)yes 

i£(yea.eq. 'y' .or.yea.eq. *Y')  then 
go  to  1030 
else  • 

75  write(*,60) 

write(*,70) 
write(*,71) 

71  format ( lx, ' (3)  Change  nozzle  data.') 

writa(*,72)  , 

72  format(lx,'  type  1,  2,  or  3  ') 

read(*,100)  1 

if (l.lt.l.and.i.gt.3)  then 

write(*,*) 'Error  iiiputJ  1  1,  2  or  3  must  be  typed* 
go  to  75 
else 

if(i.eq.l)  go  to  1000 
if(i.eq.2)  go  to  1010 
if(i.eq.3)  go  to  1020 
endif 
end  if 

c******* ***************************************** *********** 


c  * 

C  CHOOSE  THE  APPROPRIATE  GRAIN  DESIGN  * 

C  NOTE:  OPTION  2,  3,  4,  5,  6,  7  ARE  NOT  AVAILABLE.  INPUT  * 
C  YOUR  OWN  DESIGN  BY  SUBSTITUTE  THE  WRITE  AND  GOTO  STATE-  * 
C  MENTS  IN  THE  CORRESPONDING  BLOCK  * 

C  * 


C«************««***«*****tktk*******  ********************  ****** 

1030  write(*,li0) 

110  format (lx, 'Choose  one  of  the  following  grain  configurations') 
write(*,120) 

120  format ( lx, ' (1)  Star  grain') 

write(*,l30) 

130  format (lx, ' (2)  Optional  grain  ') 
write(*,l40) 

140  format(lx, ' (3)  Optional  grain  ') 
write(*,l50) 

150  formatdx, '  (4)  Optional  grain  ') 
write (*,160) 

I'O  format (lx, ' (5)  Optional  grain  •) 

write (*,170)  , 

170  format ( lx, ' (6)  Optional  grain  J) 

write (*,180) 

180  format (lx, ' (7)  Optional  grain  ') 
write (*, 190) 

190  format (lx, 'Type  1,  2,  3,  . etc.') 

read(*, *) i 

if (i.lt.l.or.i.gt.7)  then 
write(*,200) 

200  format (lx, 'Error  input! !  must  type  1,  2,  3, . ,7') 

go  to  1030 
else 

If(i.eq.l)  then 
rr5=d/2 . 0 

write(*, *) 'Run  "DESIGN  MODE"  ' 

1050  call  input2(pfac,aport,webrgd,rr5,l,tb,xx) 
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call  grains 

writs (*,*) 'Are  the  results  are  acceptable?  y-acceptable  n-not 
•acceptable  type  y  or  n. . * 
read (*,20) yes 

if (yes. eg. 'n' .or.yes.eq. 'N')  then 
go  to  1050 
else 

go  to  1055 
endif 

1055  rr5-d/2.0 

write (*,*) 'Run  the  "BURNBACK  MODE"  to  get  web  vs.  burning  area' 
call  input2 (pfac,aport,webrgd,rr5, 1, tb,xx) 
call  grains 

wrlte(*,*) 'Are  the  results  are  acceptable?  y-acceptable  n-not 
•acceptable  type  y  or  n..' 
read (*,20) yes 

if (yes. eg. 'n' .or. yes. eg. 'N' )  then 

write (*,*)' Do  you  want  to  ti.y  anot..er  star  grain?  y/n  ' 
read(*,20)yes 

if (yes. eg. 'n* .or. yes. eg. 'N')  then 
go  to  1030 
else 

go  to  1050 
endif 
else 

go  to  1060 
endif 
else 

if (i. eg. 2)  then 
writa(*,280) 

280  format (lx, 'Input  your  own  grain  design  pac)cage') 

go  to  1030 
else 

if (i. eg. 3)  then 
write(*,290) 

290  format (lx, 'Input  your  own  grain  design  package') 

go  to  1030 
else 

if (i. eg. 4)  then 
write(*,300) 

300  format ( lx, ' Input  your  own  grain  design  package') 

goto  1030 
else 

•  if (i. eg. 5)  then 
write(*,310) 

310  format (lx, 'Input  your’own  grain  design  package') 

go  to  1030  ' 

else 

if(i.eg.6)  then 
write(*,320) 

320  format (lx, 'Input  your  own  grain  design  package') 

go  to  1030 
else 

if(l.eg.7)  then 
write(*, 330) 

330  format (lx, 'Input  your  own  grain  design  package') 

go  to  1030  j  . . 

endif 
endif 
endif 
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endif 

endif 

endif 

endif 

endif 

c 

c*********************** ************************************ 

c  •  * 

c  INPUT  THE  DATA  FOR  THRUST  &  PRESSURE  VERSUS  TIME  PROFILE* 
c  CALCULATION  * 

c  * 

c*********************************************************** 

1060  write(*,*) 'The  following  processes  is  to  calculate  the  "THRUST  VS. 
*  TIME"  and  "PRESSURE  VS.  TIME"  profiles'  • 

va- (d/2.0) **2.0*3.14159*1 
rrqd"webrqd/tb 
write (*, 106 ) rrqd 

106  format (lx, 'Required  burning  rate  for  new  web,  IN/SEC..  ',f6.4) 

1070  call  rc]ctin(ath,ae,va,tamb,pa!nb,pc,atacf,rrqd,pik,etac) 
write(*,*) 'Any  change  for  the  input?  y-yes,  n-no  ' 
read (*,20) yes 

if (yes.eq. 'y' .or.yes.eq. 'Y' )  then 
go  to  1070 
else 

go  to  1080 
endif 
c 

C  * 

c  COMPUTE  THE  PRESSURE  VS.  TIME  &  THRUST  VS. TIME  PROFILE  * 

C  OF  ROCKET  MOTOR  * 

C  * 

c************************************ *********************** 

c 

1080  call  rocket 
c 

write(*,340) 

340  format (lx, 'Are  the  results  acceptable?  y-yse,  n~no . ') 

read(*,20)yes 

if (yes.eq. 'y' .or.yes.eq. 'Y')then 
go  to  9999 
endif 

write (*,350) 

350  format  (lx, 'Do  you  want  to  proceed?  y-yes,  n-no . ') 

read (*,20) yes 

if (yes.eq. 'N' .or.yes.eq. 'n')  then 

write (*,*) 'Terminated  on  the'  request  of  user  ' 
go  to  9999 
endif 

1040  write(*,360) 

360  format (lx, 'Suggest  to  change  the  following  data') 
write (*,380) 

380  format ( lx, ' (1)  Change  the  Grain  Design.') 

write(*,390) 

390  format(lx, ' (2)  Change  the  Ingredients  of  the  propellant.') 
writa(*,400) 

400  format(lx, ' (3)  Change  the  Basic  Design  Requirements.') 
write(*,*) ' (4)  Re-run  the  "ROCKET"' 
write(*,*)'  Type  1,  2  or  3  for  selection' 
read ( * , 100) i 
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if (i.ne. l.and. i.ne. 2. and. i,ne. 3)  than 

write (*,*)'  Error  input  you  have  to  type  1,  2  or  3' 
go  to  1040 
else 

if(i.eq.l)  then 
go  to  1030 
else 

if(i.eq.2)  then 
go  to  1010 
else 

if(i.eq.3)  then 
go  to  1000 
else 

if(i.eq,4)  go  to  1070 
endif 
endif 
endif 
endif 

c********* ************************************ ************** 

c  * 

c  INPUT  THE  TRAJECTORY  DATA  (FUTURE  WORK)  * 

c  * 

c* ************************************************** ******** 

c 

c  call  flyin 

c 

c* ************************************************** **^^*^^^ 
C  • 

C  CALCULATE  THE  TRAJECTORY  OF  ROCKET  KOTOR  (FUTURE  WORK)  * 

C  * 

C*********************************************************** 


call  flyit 
9999  end 


t 
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subroutine  Input 0 ( f avg , tb , 1 , d , meop , pc , r , pamb , tamb , h , f min , tpmax , 
*tpnom , tpmin , pik) 

C**********************************^' ************************ 

c  * 

c  THIS  SXJBHOOTINE  IS  FUNCTION  AS  A  BASIC  INPUT  FOR  SOLID  * 
c  PROPELLANT  ROCKET  MOTOR  DESIGN  * 

c  »  * 

c**** ************************************************  ******* 

c 

c  variable  definition: 

b  *********  (ballistic  performance)  ********* 


c  favg  :  average  thrust  (F12.6) 

c  tb  :  burning  time  of  motor  (F12.6) 

c  fmin  :  minimum  thrust  {F12.6)  • 

c  *********  (mission/vehicle)  *************** 

c  1  :  motor  length  (F12.6) 

c  d  :  motor  diameter  (F12.6) 

c  meop  :  maximum  expeted  operating  pressure  (Fl^.6) 

c  va  :  motor  volume  avaiable  for  propellant  (F12.6} 

c  tpmax,  tpnom,  tpmin 

c  pmaxtp  :  meop  *  (1-reproducibility  tolerance)  (F12.6) 

c  pc  :  norminal  chamber  pressure  (F12.6) 

c  rptol  :  reproducibility  tolerance  (F12.6) 

c  r  :  burning  rate  (F12.6) 

c  web  :  web  of  grain  (F12.6) 

c  pik  :  temperature  sensitivity  of  propellant (F12 . 6) 

c  h  ;  altitude(F12.6) 


real  favg , tb , fmin , 1 , d , meop , va , dltp , pmaxtp , pc , rptol , r , web , h 
c  open  a  file  which  can  hold  the  basic  input  data 

open ( 2 , file- ' inputo . dat ' , access- ' sequential ' , status- ' unknown ' ) 
c  input  basic  data 

write(*,*) 'All  input  data  in  decimal  format* 

write (*,*) 'Specify  the  average  thrust,  LBF .  ' 

read  (*,10)favg 
10  format(fl2.6) 

write(2, 10) favg 

write(*,*) 'Input  burning  time  of  motor,  SEC .  ' 

read(*,10)tb 
write (2, 10) tb 

write(*, *) 'Input  minimum  thrust  required,  LB .  ' 

read (*,10) fmin 

write (*,*)' Input  motor  length,  INCHES .  ' 

read (*,10)1 
write (2, 10)1 

write (*,*) 'Input  motor  diameter,  INCHES .  ' 

read(*,10)d  a 

write (2, 10) d  * 

write (*,*) 'Input  maximum  expected* • 

write(*, *) 'operating  pressure,  PSI .  ' 

read(*, 10)meop  ' 

write (2, 10) meop 

write (*,*) 'Input  Tpmax,  F .  ' 

‘  read(* , 10) tpmax 
write (2,10) tpmax 

write (*,*)' Input  Tpnom,  F .  ' 

read(*, 10) tpnom 
write (2 , 10) tpnom 

write (*,*)' Input  Tpmin,  F .  ' 

read ( * , 10) tpmin 
write (2 , 10) tpmin 
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write(2 , 10)dltp 

writa(*,*) 'Input  reproducibility  tolerance,  \ .  ' 

•  read(*,10)rptol 

write(2,10)rptol 
rptol-rptol/100.0 
pmaxtp«nneop* { 1 . 0-rptol ) 
rptol-rptol *100.0 
r-(d/4)/tb 
write(*,20)r 
write(2,10)r 

20.  format  (lx, 'Estimated  repuired  bum  rate .  ',fl2.3, 

*'  INCH/SEC •) 
write(*,30) 

c  The  following  calculation  of  ambient  temperature  and  pressure  at 
c  certain  altitude  is  directly  copied  from  KWC  Technical  Memorandum 

c  4757  subroutine  atmstd  based  on  U.5  Standard  Atmosphere, .  1976. 

30  format (lx, 'Input  the  required  design  altitude  of  rocket  motor  in  F 

*T  . 

*’) 

read(*,10)h 
write(2, 10)h 
if (h-36152.) 1,1,2 

1  tamb-S18.67-0.003559969*h 

ro»0.00237696*exp(-(0.028800339*(h*.001)+0.000132317*(h*.001) **2) ) 
go  to  7 

2  if (h-6582S.)3,3,4 

3  tamb»389.97 

ro-0. 00237696*0. 29697*exp (-0. 04771923 *(h-3 6152. ) *.001) 
go  to  7 

if (h-105511.)5,5,6 
tamb-3S4 . 142+0 . 0005442926*h 

ro-0. 00237696*0. 0718594*exp(-(0. 049130937* (h-64a25.)*. 001 
S.-0. 000034655*  ( (h-65825 . )  * .  001)  **2) ) 
go  to  7 

tamb-251. 4114+0. 00151794*h 

ro»0. 00237696*0. 0107993*exp(-(h-105511.)/21320.) 
pamb»ll. 9203*ro*tamb 
write(2,10)pamb 
write(2,10)tamb 

end  of  calculation  of  tempressure  and  temperature  at  design  alt. 

write (*,*) 'Based  on  available  requirements,  make  an  initial' 
write(*, *) 'selection  of  propellant.  Then  input  the  temperature' 

write (*,*) 'sensitivity  of  the  propellant,  %/F .  '  • 

read(*,40)pik 

write(2 ,40)plk  > 

40  fcrmat(fl2.6)  * 

pik-pik/ 100.0 

pc-pmaxtp*exp(pik* (tpnom-tpmax) ) 

pik-pik*100 . 0 

write (2, 10) pc 

close ( 2 , status- ' keep ' ) 

return 

end 
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subroutine  inputl (pc, pamb) 

c* ********************************************************** 


o'  * 
C  THE  PRIMARY  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  ESTABLISH  * 
c  The  DATA  FILE  WHICH  CAN  BE  ACCESSED  BY  MICROPEP  SUB-  * 
c  TINE  * 

c  * 


c* ************************************** ******************** 

c. 

common  /mlcrp/  a(12,12) ,Xr(20) ,amat(10,12) ,jat(12) , in, Is, fie (10, 6) 
l,le(10,6) ,alp(12) ,w27,n,dh(10) ,rho^lO) ,wate(10) ,wl(6) ,w43,lg,np, 
2vnt(201) , W47, name, ser,floor,ltag (100) ,wlng(10) 
common/ chara/block , j  e , aspec , specie 

character's  specie (200)  , 

character*30  block ( 10 ), of lie 
character*2  je(10,6) ,aspec(12) 
common/moon/tstest , te , Irun , lounlt , lend 
''haracter*30  blok(lO) 
character*l0  user 
dimension  jle(10,6) 

Integer  nlc,niu, irun, alt, kr 
real  denexp 

open  the  file  as  formatted  sequential  as  unit  3 
open (3 , file- ' lnput.dat ' , access*' sequential ' , status- ' unknown ' ) 
c  Input  the  case  name  from  the  keyboard 

write(*,*) '  Start  to  Input  the  selected  propellant  data  ' 

write (*, *)  '  ' 

write(*,l) 

1  format (lx, 'Specify  the  case  name  (format  AlO) .  ') 

read(*, ' (a) ' ) user 
write (3, 10) user 
10  format (lalO) 

irun-2 

c  Input  the  number  of  Ingredients,  supplied  ingredients,  and  runs 
write (*,•)' Input  the  number  of  ingredients?  (format  15)....' 
read(*,*)nic 

write (*,*) 'Input  the  number  of  user  suppplied  ingredients  (15)...' 
read(*, *) niu 
wr ite (3,25)nic,niu, irun 
20  format(2i5) 

25  f=r=at(3i5) 

c  input  density  exponent  which  is  used  in  conjunction  with  control 
c  option  3. 

write  t*i*)'Input  the  density  exponent;' 

write(*,*) 'For  tactical  missiles  or  first  stages  use  1.0' 

write (*,*) 'Use  0.7  for  second  stagps  and  0.2  for  third  stages' 

write(*,*)  '  (format  FIO.O) . . . .  ' 

read (*, 4 0 ) denexp 
write (3,40) denexp 
40  format(fl0.5) 

c 

c  input  control  option  store  it  in  kr(20) 
c  print  out  the  options  on  the  screen 

writa(*, ' (a) ' ) 'There  are  five  options  listed  below  for  your' 

write(* ,*)' selection.  Type  1,  2,  3,  4,  5 - ' 

write(*, *) 'Normal  input  for  design  is  3' 

write(*, ' (a) ' )  '  1.  For  chamber  and  exit  calculations.' 

write(*, ' (a) *) •  2.  For  chamber  only  calulations. ' 

write{*, ' (a) ' ) '  3.  For  boost  performance  and  nozzle  design.' 
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write (* , ' (a) ' ) '  4.  For 

write(*, ' (a) ') '  and 

write (*, ' (a) ') '  in 

write (*, ' (a) ’) •  5.  0th 

write(*, ' (a) ' ) '  alt 

write(*,*) 'Use  format  II 
read(*,45)alt 
format ( il) 
if  (alt. It. S)  then 
kr(2)  -  0 
lcr(4)  -  0 
kr(5)  -  0 
kr(6)  -  0 
kr(8)  -  0 


4.  For  equilibrium  calculations  at  specified  T' 
and  P.  Use  T  in  Kelvin  at  Pe  location  and  P' 
in  psia  at  the  normal  Pc  location. ' 

5.  Other  options  (  Check  README  for  special* 
alternatives) . • 


if  (alt.eq.l)  then 
kr(l)  -  0 
kr(3)  -  0 
kr(7)  «  0 

else 

if  (alt.eq.2)  then 
kr{l)  -  1 
kr(3)  -  0 
kr(7)  -  0 

else 

if  (alt.eq.3)  th 


kr(l) 

lcr(3) 

kr(7) 


lcr{l) 

kr(3) 

)cr(7) 


endif 
end  if 
endif 

write (3, 50) (kr(i)  ,i 
format (8il) 


else 

write (*, ' (a) ' ) '  Input  your  requirements  for  calculation--' 


read(*,60) (kr(i) ,i  -  1,15) 
write(3,60) (kr(i) ,i  -  1,15) 
format  (15il) 


endif 

input  the  serial  number  of  ingredients  found  in  JANNAF.DAT  table 


of  number 


Ingredient (15) . . . 


write(*,70)  nic 

format (lx, 'Input  ', i2 ,'■ ingredients  needed  for  calculation') 
do  100  i  «  l,nic  •  . 

write(*,80) i 

format(lx, 'Serial  #  of  number  ',i5,'  ingredient (15) .. .  ' 

read (* , 90) itag ( i) 

format(i5) 

3  continue 

write  (3, 120) (itag(i) ,  1  -  l,nic) 

3  format(10i5) 

Input  user  supplied  ingredients 
if  (niu.gt.O)  then 
do  170  i  -  l,niu 

write(*, ' (a) ') 'Specify  ingredient  name  use  caps  lock  (A30) 
read(*, 130) blok (i) 


l,nic) 
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write(*r' (a) ') 'Specify  composition  of  ingredient  max  »  6' 
write(*, ' (a) •) •  use  (6(13, A2))' 
write(*, ' (a) ') •  ###**###**###**###**###**###**• 
read(*,140) (jie(i,l) , je(i,l) ,1-1, 6) 

wrlte(*, ' (a) ' ) '  Specify  heat  of  formation,  and  density  (f5.0,f6 
*.0)  • 

read(*, 15Q)dh(i) ,rho(i) 

write(3,160)blok(i) , (jie(i,l) ,je(i,l) ,1-1,6) ,dh(i) ,rho(i) 
format (a30) 
format (6 ( i3, a2 ) ) 
format (f 5 . 0 , f 6. 0) 
format (a30, 6 (i3,a2) ,fS.0,f 6.0) 
continue  ’ 

else 

write  (*,'(a)')'  No  user  supplied  ingredients  ' 
endif 

if  (irun.lt. 2)  then 
irun  -  2 
endif 
status-0 

do  220  i  -  l,irun 

if (status. eg. 0)  then 
wl(5)-2S0 

write  (*,'(a)')'  Specify  the  desired  chamber  pressure  (F12.4)  ’ 
read  (*,180)wl(5) 
wl (6) -pamb 

write  (*,'(a)')'  Specify  the  exhaust  pressure  (F12.4) ' 

read  (*,190)wl(6) 

ni-nic+niu 

do  175  j  «l,ni 

write(*,172) itag(j) 

format (lx, 'Input  the  weight  %  of  ingredient i5, ' (FIO. 4) ..." ) 

read(*, 190) wing( j ) 

continue 

format(fl0.4) 

format(fl0.4) 

write(3,210)wl(5) ,wl(6) , (wing(j) , j-1,10) 

format(12fl0.4) 

status-1 

else 

wl(5)“pc 

write(3,210)wl(5)  ,wl(6) ,  (wing(j)  ,  j-1',10) 
endif ■ 
continue 

close (3)  ? 

return  *  ^ 

end 
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subroutine  pf e { f avg , tb , 1 , d , pc , cf c , ath , etacf , m , vp , vl , dpe , aport , 

+pf  ac , isp , c  f th , webrqd , wf , tc , r , aa , etac , xx ) 
g*********************************************************** 
c  * 

c  THE  MAIN  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  CALCULATE  THE  * 

C  PERFORMANCE  LOSS  AND  THE  DATA  FOR  GRAIN  DESIGN  * 

c  * 

c*********** **************************************** ******** 

c  variable  definitions: 

c  favg:  Average  thrust  from  basic  input 

c  cfth:  thrust  coefficient  from  micropep 

c  cfc  :  corrected  cf 

c  pc  ;  chamber  pressure  from  basic  input 

c  d  :  motor  diameter  from  basic  input  > 

c  1  :  motor  length  from  basic  input 

c  athth:  theoretical  throat  area 

c  ath  :  throat  area 

c  atht  :  tempory  result  of  throat 

c  dthth:  theoretical  throat  diameter 

c  etacf:  thrust  coefficient  efficiency 

c  dif  :  for  convergence  test  of  throat  size 

c  va  :  vol  avaiable  for  propellant 

c  tb  :  burn  time 

c  vp  :  propellant  vol. 

c  vl  :  volumetric  loading  factor 

c  aport:  port  area 

c  dpe  :  equivalent  circular  port  diameter 

c  pwet  :  wetted  perimeter 

c  j  :  throat-to-port  area  ratio 

c  rhop  :  propellant  mass  density 

c  mw  :  molecular  weight  of  propellant 

c  mf  :  mole  fraction  of  condensed  phase 

c  webrqd:  required  web 

c  wf  :  web  fraction 

a  pfac  :  perimeter  factor 

c  character* 80  fnama 
c 

real  favg , cfth , cf c , pc , d , 1 , athth , ath , dthth , etacf , dif , tb , vp , vl , va , j , 
+aport , dpe , rhop , a , nn , rstd , rpc , web , webrqd , wf , m , pwet , pf ac , isp , ispm , at 
+ht,mw,tc,ispth,ld,mf 
character*!  yes 

c  The  purpose  of  this  common  instruction  is  to  transfer  w43 (propellant 
c  density)  from  MICROPEP;  SUBROUTINE  PUTIN(LE) 

common  /micrp/  a(12,12) ,kr(20) ,amat(10,12) , jat(12) ,in,is, 
lfie(10,6) ,ie(10,6) ,alp(12) ,w27,n,dh(10) ,rho(10) , 

2wate(10) ,wl(6) ,w43,ig,np,vnt(20li,w47,name,ser 
3,floor,itag(100) ,wing{10)  » 

C  The  purpose  of  this  common  instruction  is  to  transfer  the  CHAMBER  PRE- 
C  SSURE  S  MOLECULAR  WEIGHT  from  MICROPEP  SUBROUTINE  OUT. 
common  /tcf/  tfc,rmwtc 

C  Transfer  SPECIFIC  IMPULSE,  CHARACTERISTISC  VEL. ,  CODFF.  OF  THRUST  & 

C  OPTIMUM  EXPANSION  RATIO  from  MICROPEP  SUBROUTIN  DESING. 

common  /socf/  spi (2) ,cf (2) ,oex(2) 
c  calculate  the  theoretical  throat  area 
pi-3.14159 

c  transfer  cfth  from  "MICROPEP" 
cfth-cf(2) 
c  initialize  cfc 
cfc-cfth 

athth=favg/ (cfth*pc) 
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c  initialize  ath-athth 
ath-athth 

o  calculate  theoretical  throat  diameter 
dthth»sqrt (ath/pi) *2 
write ( * , 88  8 ) cf th , athth , dthth 

888  fonnat(lx, 'Cfth-' ,f5.3, '  Athth-' ,f 6.3, '  Dthth- ' , C5 . 3 ) 
c  calculate  the  etacf 

write(*,*) 'Input  the  nozzle  half  angle' 

write(*,*) 'alpha,  DEGREES  (F5.2) .  ' 

read(*, *) alpha 
pi-3.14159 

alpha=(alpha/iao.O) *pi 

writa(*, *) 'Input  the  nozzle  exit  half  angle' 

write(*,  *) 'thetaex,  DEGREES  (F5.2) .  '  , 

read ( * , * ) thetaex 
thetaex=(thetaex/180.0) *pi 

write (*,*)' Input  radial  erosion  rate  of  the  throat  in  IN/SEC 
write {*,*)' assumed  to  sLart  at  time  zero  (FORMAT  O.XXXX)....' 
read(*, *)erosr 

write (*,*) 'Input  the  mole  fraction  of  condensed' 

write(*,  *■  'phase  in  (MOLE/100  g) . . .  ' 

read(*, *)mf 

write (*,*)' Input  the  material  of  nozzle:' 

write (*,*)' 1-steel  nozzle,  0-ordinary  nozzle .  ' 

read  (*,*)nozle 
atht-ath 
do  9  i-1,50 

C  calculate  the  loss  mechanism 

call  loscf (etacf ;tb,ath, pc, alpha, thetaex, ero3r,mf,nozle) 
c  correct  the  trust  coefficient 
cfc-cfth*etacf 

c  calculate  the  corrected  throat  area 
ath-favg/ (cfc*pc) 
write (*, 8888 ) etacf ,cfc,ath 

8888  format (lx, 'Etacf-' ,fl2. 3, '  Cfc-' , fl2. 3 , '  Ath-',fl2.3) 

c  calculate  the  convergence  athth  and  ath 
dif-abs ( (atht-ath)/atht) *100.0 
if  (dif.gt.1.0)  then 
atht-ath 
else 

go  to  91 
endif 

9  continue  , 

91  va-pi* (d/2) **2*1 

etac-0.93 
vol-pi* (d/4) **2*1 

c  Density  and  molecular  weight  com  from  "MICROPEP" 
rhop»w43 
mw-rr.wtc 
tc-tfc 
ispth-spi (2) 

write (* , 98) rhop,mw, tc, ispth 

98  formate  Rhop- ' , f 12 . 3 , '  MW-',fl2.3,'  Tc-',fl2.1,'  Ispth- ' , f 12 . 2 

*) 

99  i3p=ispth*etaof *etac 
m=favg/isp 

vp=f avg*tb/ ( isp*rhop) 
vl-vp/va 

write(*,87)VP,VA,VL 

87  formate  VP=',F12.1,'  VA-',F12.1,'  VL-',F12.3) 
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dpe-sqrt ( 1. 0-vl) 
aport= (pi/4 ) *d**2* ( l-vl ) 

tres- ( aport*l) *pc*mw/ (m* 12 *1546* (tc+459 . 6) ) 


tres“tres*1000 . 0 
write (*,88) tres 

88  formate  RESIDENCE  TIME  IN  MSEC.,. .  •,fl2.3) 

write (*, 89 )etac 

89  format (*  COMBUSTION  EFFICIENCY  (etac) .  ',fl2.3) 


write(*,  *) 'Change  etac?  y-yes,  n-no .  ' 

read (*,20) yes 

if (yes.ge. 'y' .or. yes. eg. 'y')  then 

write  (*,*)' Input  etac .  ' 

read (*,30) etac 

vol“aport*l  ‘ 

go  to  99 
endif 

c  calculate  the  nozzle  exit  area 
ae=ath*oex(2) 

c  calculate  the  required  perimeter  factor 
write (*, 19) pc 

19  format ( lx, ' Input  the  C*  with  Penom*' , flO. 1, '  (F10.4) .  ') 

read(*, *) cstar 

xx= (pc*ath*32 . 17) / (cstar*rhop*l*3 . 14159*dpe) 

18  write (*,*)' Input  the  estimated  required  web  (INCHES)  (F10.4)' 

write (*,*)' suggest  initial  guess  D/4 .  * 

read ( * , * ) webrqd 
r=webrqd/ tb 
pfac*xx/r 
write(*,21) r 

21  format (lx, 'The  required  burn  rate  »',f6.4) 
write ( * , 22) pfac 

22  format (lx, 'The  required  perimeter  factor*’ , f 6. 3) 
write (*,25) XX 

25  format(lx, 'pfac*' ,f5.4, '/r') 

write(*,24)tb 

24  format(lx, 'tb  »',f7.2) 

wr ite ( * , *) ' Are  these  acceptable?  y-yes  n-no .  ' 

read ( * , 20} yes 

if (yes.ge. 'y' .or.yes.eq. 'Y' )  then 
pwet*pf ac*pi*dpe 
wf=webrqd/ (d/2 ) 
j=ath/aport 
ld=l/d 
go  to  23 
else 

go  to  18  3 

endif  i 

20  format(al) 

30  format(fl2.6) 

23  write (*, 40) cfc 


40  format  (lx,  'THRUST  COEFFICIENT . ',F12.3) 

write ( * , 50) etacf 

50  format  (lx,  'ETACF . ',F12.3) 

write ( * , 60) isp 

60  formatdx,  'SPECIFIC  IMPULSE,  LBFSEC/LBM . ',fl2.3) 

write(*, 70)m 

70  format  (lx,  'MASS  FLOWRATE,  LBM/SEC . ',fl2.3) 

write  f  * , 80) vp 

80  format  (lx,  '  PROPELLAJIT  VOLUME,  CUIN . ',fl2.1) 

write ( *, 100) ath 
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format ( 
write (* 
format ( 
write (* 
format  ( 
write (* 
format  { 
write { * 
format  ( 
write (* 
format ( 
write ( * 
format  ( 
write ( * 
format  ( 
write  ( * 
format ( 
write (* 
format ( 
write ( * 
formate 
write  (* 
format ( 
return 
end 


lx, 'NOZZLE  THROAT  AREA,  SQ  IN . 

,105)ae 

lx,  'NOZZLE  EXIT  AREA,  SQ  IN . '  , 

,110)dpe 

lx, 'EQUIVALENT  PORT  DIAMETER,  IN . ' , 

,120)aport 

lx,  'PORT  AREA,  SQ  IN  . ', 

,130)pfac 

lx,  'PERIMETER  FACTOR . ' , 

, 140) webrqd 

lx,  'REQUIRED  WEB,  IN . ', 

,150)wf 

lx,  'REQUIRED  WEB  FRACTION . '  , 

,90)vl 

lx,  'VOLUMETRIC  LOADING . . . '  , 

,190)ld 

lx,  'LENGTH  TO  DIA!-IETER  RATIO . ', 

, 170) tres 

lx, 'RESIDENCE  TIME,  MSEC . ’, 

,180)  j 

lx, 'TiiROAT  TO  PORT  AREA  RATIO  (J  FACTOR) . . . ' , 
,200)r 

lx,  'REQUIRED  BURNING  RATE,  IN/SEC . ', 


fl2.3) 

fl2.3) 

fl2.3) 

fl2.3) 

fl2.3) 

fl2.3) 

’fl2.3) 

fl2.3) 

fl2.3) 

fl2.3) 

fl2.3) 

fl2.3) 
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subroutine  loscf ( etacf , tb , ath , pc , alpha , thetaex , erosr , 

*fflf ,nozle) 

c* *****************«**««*««*******•• ************************ 

c  * 

c  THIS  SUBROUTINE  CALCULATE  THE  NOZZLE  LOSS  COEFFICIENT  * 

c  CTHRUST  COEFFICIENT)  USING  THE  METHODS  PROPOSED  BV  AGARD  * 
o  PROPULSION  AND  ENERGETICS  PANEL  FOR  METALLIZED  PRO-  * 

c  PELLANTS  (Reference:  Performance  of  Rocket  with  Metal-  * 
c  lized  propellants,  AGARD-AR-230)  ♦ 

c  * 

C*********‘************************* ************************* 

c  Defination  of  varibles  in  this  subroutine 
c  etacf:  nozzle  loss  coefficient 

c  div  :  percent  divergence  loss  • 

c  kin  ;  percent  kinetics  loss 

c  bl  :  percent  boundary  layer  loss 

c  tp  :  percent  two-phase  loss 

c  sub  :  pcrsent  submergence  loss 

c  eros  :  percent  nozzle  erosion  loss 

c  astar:  nozzle  entrance/  nozzle  throat  area 

c  s  :  submergence  length/  length  of  internal  motor 

c  mf  :  mole  fraction  of  condensed  phase  in  mole/lQQgrm  • 

c  cl,c2,c3,c4,c6  : 

c  constant  in  emperical  relation 

c  dp  :  diameter  of  particle 

c  pc  :  chamber  pressure 

o  ispf  :  isp  of  frozen  equilibrium 

c  isps  :  Isp  of  shifting  equilibrium 

c  nozle:  steel  nozzle  =•  1,  ordinary  nozzle  »  0 

c  tb  :  burning  time  of  motor 

c  dt  :  throat  diameter 

c  expn  :  nozzle  expansion  ratio 

c  ispi  :  Isp  with  ideal  (at  design  point)  nozzle  expansion  ratio 

c  ispm  :  Isp  with  mean  nozzle  expansion  ratio 


integer  expnu,expnl 

real  etacf , div, kin, bl , tp, sub, eras, astar, s,mf, dp, pc, isps, ispm, ispf, 
+dt , ispu , ispl , erosr 
common  /scratc/  plot(5,100) 
common  /socf/  spi (2) , cf (2 ) ,oex(2) 
c  calculate  the  %  of  nozzle  divergence  loss 

div  ■  50  *  (  1  -  cos  ((alpha  +  thetaex)  /  2  )) 
c  ispf  is  the  frozen  and  isps  is  the  shifting  specific  impulse  from 
c  micropep  run  'at  nominal  pc 

c  Transfer  specific  impulse  data  from  MICROPEP.' 
ispf»spi(l)  i 

isps“spi(2) 

c  Calculate  the  kinetic  loss 

kin  -  33.3  *  (1  -  ispf/ isps) 
c  calculate  %  boundary  layer  loss 

c  Steel  nozzle  (1),  Ordinary  nozzle  (0).  Enter  0  or 
if  (nozle. eg. 1)  then 


cl 

c2 

else 

cl 

c2 

endif 


0.00506 

0.0 

0.00365 

0.000937 
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c  expansion  ratio  of  the  nozzle  for  shifting  equilibrium  from 

c  micropep  with  exit  pressure»ambient  pressure  at  design  altitude 

c  read  (*,70)expn 
expn»oex(2) 
pi=3. 14159 
dt»sqrt(ath/pi) *2 

bl  ■(cl*pc**0.8/ (dt*0.2) ) *(l+2*exp(-c2*pc**0.8*tb/(dt*0.2) ) ) *(1+0. 
$016* (expn-9. 0) ) 

c  calculate  the  two-phase  flow  loss 

c  the  mole  fraction  of-  condensed  phase  is  come  from  "MICROPEP" 
c  mean  particle  diameter  calculation  by  using  empirical  relation 
dp  ■  3.39*dt**0.4692 
if  (mf.lt. 0.09)  then 

mf  «  0.09  ' 

endif 

if (dt.lt. 1.0)  then 
c3=9.0 
C5-1.0 
C6-1.0 

else 

if(dt.le.2)  then 
C3-9.0 
C5-1.0 
C6=0.8 

else 

if (dp. It. 4)  then 
C3-13.4 
c5»0.B 
C6=0.8 

else 

if(dp.le.S)  then 
C3-10.2 
C5«0.8 
C6-0.4 

else 

C3-7.58 

C5-0.8 

c6=0.33 

endif 

endif 

endif 

endif 

C4-0.5 

tp»c3*mf**c4*dp**c5/ (pc**0. 15*expn**0.08*dt**c6) 
c  calculate  the  submergence  loss (NOT  USED  IN  PRELIMINARY  WORK  U) 
c  the  nozzle  is  assumed  to  be  of  external  design  with  no  submergence 
c  loss.  s»length  of  submergence  /  length" of internal  motor, 
c  astar»‘(pi*d**2/4)/ath 

c  sub»0. 0684* (pc*expn/astar) **0 . 8*s**0. 4/dt**0. 2 

c  default  to  0.0 

sub^O.O 

c  calculate  the  nozzle  erosion  loss 

c  calculate  the  ispm  (specific  impulse  at  mean  expansion  ratio) 
df-dt+erosr*tb*2 
dbar={dt+df)/2 
expnm»'(dt**2/dbar**2)  *expn 
expnl»int ( expnm) 

c  Transfer  expansion  ratio  of  nozzle  and  specific  impulse  from  MICROPEP. 
ispl«plot ( 4 , expnl) 
expnu=int (expnm) +1 


■  ^  ^  ^  ^  «  -  •  _  •  _  ^  ^  .  •  *.  •.  M  ~  m.  *  a  *  «■  ~  ^  j 
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ispu“plot ( 4 , expnu) 

ispm= (expnm-expnl)/ (expnu-expnl) * (ispu-ispl) +ispl 
.c  erosion  loss  calculation 

eros-(l-ispni/isps)  *100 
write (*, 100) div 

100  format  (lx, 'PERCENT  DIVERGENCE  LOSS .  ',f6.3) 

write (*, 110) tp 

110  format  (lx, 'PERCENT  TWO  PHASE  FLOW  LOSS .  ',f6.3) 

write (*, 120) bl 

120  format  (lx, 'PERCENT  BOUNDARY  LAYER  LOSS .  *,f6.3) 

write (*,130) kin 

130  format  (lx, 'PERCENT  KINETIC  LOSS . ,',f6.3) 

write(*, 140) sub 

140  format  (lx, 'PERCENT  SUBMERGENCE  LOSS .  •,f6.3) 

wr ite ( * , 1 5  0 ) eros 

150  format  (lx, 'PERCENT  EROSION  LOSS .  ',ff  3) 

etacf=l-  (div+tp+bl+kin+sub-t-eros)  /lOO .  0 
write ( * , 160) etacf 

160  format  (lx, 'THRUST  COEFFICIENT  EFFICIENCY .  ',f6.3) 

return 
•  end 


i 
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subroutine  input2 (pf ac , aport , web , rr5 , Igrain , tb , xx) 

c******************************** *************************** 

c.  * 

c  THE  PRIMARY  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  ES-  * 

c  TABLISH  THE  DATA  FILE  FOR  THE  CSD  "GRAINS"  PROGRAM  * 

c  * 

c************ *********************************** ************ 

character*SO  tapes 
character*!  yes 
real  Igraln, lspoke,kspoke,RR5 
integer  nmax,nspoke,ibb,ipmt 

write(*,*) 'The  following  will  be  the  data  input  for  ' 
write(*,*) 'spoked  star  grain  designs  with  radial  sides  ' 

open(4, file-'grainin.dat* , status- ' unknown ' )  , 

write(*, *) 'Run  DESIGN-1  mode  or  BURNBACK-0?  Type  1  or  0... 
read(*,*)i 
if (i. eg. 1) then 
11  write {*, 40) web 

14  r-web/tb 

pfac«xx/r 
wrlte{*,13)r 

13  format  (lx, 'Required  burning  rate,  IN/SEC .  ', 

write (*, 50) pfac 

write(*,*) 'Do  you  want  to  change  the  web?...y/n . 

read(*,5)yes 

if (yes.eq. ' Y' .or.yes.eq. 'y* )then 

write(*,*) 'Input  the  required  web,  IN . 

read (*,*) web 
go  to  14 
endif 

write (*,*) 'Max  number  of  spokes?. (13) . 

read ( * , * ) nmax 

write (*,*) 'Initial  number  of  spokes?. (13) . 

read(*,*)nspoka 

ibb-1 

iprnt-0 

wr ite ( 4 , 10 ) nmax , nspoke , ibb , iprnt 
write (*,30) aport 

30  format (lx, 'Required  port  area,  SQ  IN . 

*.3) 

write (*,40) web 

40  format (lx, 'Initial  web  thickness  in  IN . . . 

*.3) 

write(*, 50) pfac 

50  format (lx, 'Required  perimeter  factor  . 

*•3)  , 

write(*, *) 'Input  Kspoke  -  Abstract  multiplier  on  ' 
write(*,*) 'PHI2  -  if  KSP0KE=0.0  PHI2-0.0  where:' 
write(*,*) 'PHI2=ATAN(KSPOKE*WEB/RR3/COS(PHI3) )  (0. 0-1.0) 

read ( * , *) kspoke 
write ( * , 60) rr5 

60  format (lx, 'Grain  outer  radius, INCHES . 

*.3) 

write (4,20) aport, web, pfac, kspoke, rr5 

write(*, *) 'Input  initial  grain  inner  radius,  INCHES  (F7.3) 
read(*,*)ril 

write (*,*) 'Input  initial  value  of  BRl  -  spoke' 

write (*,*) 'top  corner  radius  in  IN  (Default-0.25) . 

read(*, *) sril 

write (*,*)' Input  initial  value  of  SR2  -  spoke' 
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write(*, *) 'base  coner  radius  in  IN  (Default«0 . 25) .  ' 

read(* , ♦) sri2 

write (4, 20) ril,sril,sri2 

write(*,*) 'Input  increment  on  SRI  in  IN  -  cannot  be  0.0...  ' 
read(*, *)dell 

write (*,*) 'Input  increment  on  RRl  in  IN  -  negative  value..  ' 
read (* , *) del2 

write(*, *) 'Input  increment  on  NSPOKE  Default  -  1.0 .  ' 

read(*, *)del3 

write (*,*) 'Input  increment  on  SR2  in  IN  -  cannot  be  0.0...  ' 
read(*, *)del4 

write(4,20)dell,del2,del3,del4 

write (*,*) 'Input  tolerance  on  grain  angle  sum  total  -  %' 

writa{*,*) 'of  360.  (0.0-1. 0) . . .  ' 

read(*, *)cirtol 

write (*,*) 'Input  tolerance  on  port  area  -  %  of  Aport  ' 

write(*,*)  '  (0.0-1. 0) .  ’ 

read'*  ♦)pttol 

write (*,*) 'Input  tolerance  on  perimeter  -  %  of  perimeter..  ' 

writa(*,*)  '  (0.0-1. 0) .  ' 

read ( * , * ) pertdl 

write (4,20) cirtol , pttol , pertol 

write(*,*) 'Input  Max  of  SRI  in  IN  .(F10.5) .  ' 

read(*, *) f ini 

write(*,*) 'Input  Min  of  RRl  in  IN  .(F10.5) .  • 

read(*, *) f in2 

write (*,*) 'Input  Max  of  SR2  in  IN  .(F10.5) .  ' 

read (*, *) fin4 

write (*,*)' Input  Max  of  web  in  IN  .(F10.5) .  ' 

read(*, *) f inweb 

write (*,*) 'Input  increment  for  web  in  IN.  (F10.5) .  ' 

read(*,*)webstp 

write (4,20) f ini , f in2 , f in4 , f inweb, webstp 
else 
nmax-0 

write (*,*) 'Number  of  spoices?. . .  (13) . . .  ' 

read ( * , * ) nspoke 
ibb— 1 

write(*, *) 'Do  you  want  optional  output?  y-yes,  n-no .  ' 

read ( * , 5) yes 

if (yes. eg. ' Y ' . or. yes. eg. 'y ' ) then 
iprnt=-l 
else 

iprnt=0 

endif 

write (4,10) nmax , nspoke , ibb , iprnt 
write(*,70)  i 


70  format  (lx, 'Input  initial  port  area,  SQ  IN..(F10,5) .  ') 

read (*,*) aport  ■ 
write ( * , 80) 

80  format  (lx,  '  Input  initial  web  thickness,  IN.(F10.5) .  ') 

read(*, *)web 
write (* , 90) 

90  format  (lx,  '  Input  the  perimeter  factor. .  (FIO .  5) .  ') 


read(*,*)pfac 
kspoke=0. 0 
write(*,60)rr5 

write (4 , 20) aport, web, pfac,kspoke, rrS 
write(*,*) 'Input  initial  RRl,  INCHES . (FIO . 5) 
read (* , *) ril 


APPENDIX 


write(*,*) 'Input  initial  value  of  SRI  (Default-0.25  IN.).... 
read(*, *) sril 

writa(*, *) 'Input  initial  value  of  SR2  (Default-0.25  IN.) . 

read(*,*)sri2 

write(4,20)ril,sril,sri2 

dell-0.0 

del2-0.0 

del3-0.0 

del4-0.0 

write(4,20)dell,del2,del3,del4 
cirtol-0. 0 
pttol=0. 0 

pertol-0.0  I 

write (4,20) cirtol , pttol , pertol 

f inl-0. 0 

fin2-0.0 

fin4-0.0 

f inweb-0.0 

webstp=0.0 

write (4, 20) finl, fln2 , f in4 , f inweb, webstp 

write (*,*) 'Input  grain  length,  INCHES  (F10.5) .  ' 

read(*, *) Igrain 

write(*,*) 'Input  initial  PHI2,  DEGREES  (F10.5) .  ' 

read(*, *) dphi2 

write(*,*) 'Input  initial  PHI4,  DEGREES  (F10.5) .  ' 

read  (*,*)dphi4 

write (4, 20) Igrain, dphi2,dphi4 

write(*,*) 'Input  #  of  steps  in  decimal  format' 

write (*,*)  'Input  #  of  burnback  steps  for  SR1>0 .  ' 

read(*,*)stepl 

write(*, *) 'Input  #  of  burnback  steps  for  PHI2>0  6  SRl-0....  ’ 
read(*, *) step2 

write(*,*) 'Input  #  of  burnback  steps  for  LSPOKE-0  &  PHI2>0.  ' 
read(*,*)stpp3 

write(*,*) 'Input  #  of  burnback  steps  for  ' 

write(*,  *) 'triangular  spoke  burnback . . .  ' 

read(* , *) step4 

write(*,*) 'Input  i  of  burnback  steps  for  ' 

write(*,  *) 'remaining  fuel  burnout . . .  ' 

read(*,*) Steps 

wr ite ( 4 , 2 0 ) stepl , step2 , step3 , step4 , Steps 
endif  . 

5  format (al)  ^ 

10  format (4i3) 

20  format  (5f  10 . 5)  ■, 

close (4, status-' keep ' )  a 

return 
end 
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